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2.0  AQUANAUTICS  POWER  SOURCE 


2.1  Concept 

As  part  of  DARPA's  "Artificial  Gill"  program,  Aquanautics  is  developing  an  underwater 
power  source  that  is  at  least  2  to  3  times  more  efficient  than  lithium  thionyl  chloride.  It 
uses  an  ambient  pressure  solid  polymer  electrolyte  (SPE)  hydrogen/oxygen  fuel  cell  with 
an  artificial  gill  that  extracts  dissolved  oxygen  from  seawater  and  an  Alwatt  aluminum 
corrosion  cell  that  generates  the  hydrogen.  The  'particulaHnnovation  in  this  concept  is  a 
direct  carrier  feed  fuel  cell  (CFFC)  process  which  allows  the  fuel  cell  to  operate  at  ambient 
pressure  without  a  costly  deep  water  pressure  vessel.  '  >  ^  c  '  :  / » ,  .-.u 


Figure  2.1  shows  the  schematic  diagram  for  this  system: 


|  In  this  scheme  three  basic  units  are  required  which  are  discussed  in  detail  below: 

1.  Artificial  Gill 

2.  Hybrid  Fuel  Cell 

3.  Alwatt  Hydrogen  Generator 

I 
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THE  ARTIFICIAL  GILL 

The  objective  of  the  artificial  gill  process  is  to  extract  and  concentrate  the  oxygen  that  is 
dissolved  in  the  seawater  rather  than  split  the  water  molecule  into  its  oxygen  and  hydrogen 
components  (electrolysis).  The  later  process  requires  more  energy  than  can  be  recovered 
from  burning  the  oxygen  in  a  combustion  process.  The  artificial  gill  requires  only  about 
10%  of  the  power  that  can  be  generated  from  the  extracted  oxygen  in  a  combustion 
process.  This  factor  makes  a  gill  the  the  only  feasible  alternative  for  on-site  oxygen 
production  involving  a  power  system. 

The  driving  force  for  oxygen  diffusion  is  the  difference  in  oxygen  concentration  between 
the  seawater  and  carrier.  A  membrane  is  used  to  keep  these  fluids  separated.  When  the 
oxygen  gas  diffuses  across  the  membrane  from  the  seawater  into  the  carrier  fluid,  it  binds 
to  an  oxygen  carrying  molecule.  In  this  bound  state,  the  carrier  fluid  does  not  have  any 
oxygen  partial  pressure.  Therefore,  the  driving  force  for  the  oxygen  diffusion  is  at  its 
maximum,  equal  to  the  partial  pressure  of  the  oxygen  in  the  seawater. 

Worst  case  design  conditions  are  used  for  the  gill  which  assume  oxygen  dissolved  in 
seawater  at  a  concentration  of  approximately  1  to  3  ml  per  liter  of  seawater  (see  Appendix 
G  for  measured  values  in  all  the  oceans). 

Proper  packaging  of  the  membrane  is  essential  to  provide: 

i)  low  volume  per  unit  quantity  of  oxygen  transferred 

ii)  low  seawater  drag  losses 

iii)  low  carrier  pumping  cost 

Fortunately,  the  chemical  industry  has  developed  a  wide  range  of  membrane  polymers  that 
will  work  for  underwater  systems.  Analysis  by  Aquanautics  has  concluded  that  a 
rectangular  gill  package  with  a  hollow  fiber  membrane  will  be  the  best  design  choice  with 
carrier  flowing  inside  the  hollow  fiber  and  seawater  flowing  outside.  The  details  of  this 
analysis  can  be  found  in  Appendix  B. 

THE  FUEL  CELL 

Fuel  cells  are  very  efficient  at  converting  the  chemical  energy  stored  as  hydrogen  and 
oxygen  into  electrical  energy.  The  problem  with  using  fuel  cells  at  great  water  depth  is 
finding  an  efficient  means  of  storing  the  gaseous  reactants  and  releasing  them  into  high 
ambient  pressures.  A  solution  is  to  place  the  fuel  ceil  in  a  pressure  vessel.  But  at  6,000 
meters  such  an  enclosure  would  be  very  heavy  and  expensive.  This  problem  has  limited 
fuel  cells  from  serious  consideration  for  ocean  bottom  applications. 

The  use  of  an  artificial  gill  poses  a  different  kind  of  problem.  Oxygen  dissolved  in 
seawater  can  exist  in  equilibrium  with  the  maximum  partial  pressure  of  0.2  atmospheres  or 
3  psi.  Extracting  oxygen  in  its  gaseous  form  and  compressing  it  to  above  ambient  pressure 
would  involve  a  large  expenditure  of  power  which  defeats  the  purpose  of  extraction.  In 
order  to  enable  efficient  operation  at  ambient  pressure,  it  is  desirable  to  use  the  oxygen 
while  it  is  still  in  its  bound  state  and  skip  the  intermediate  step  of  converting  it  into  a  gas. 

The  Aquanautics'  gill  process  employs  oxygen  carrying  molecules  which  can  be 
electrochemically  modulated  to  release  the  bound  oxygen.  Since  a  fuel  cell  is  also  an 
electrochemical  system,  Aquanautics  has  shown  that  the  oxygen  can  be  reduced  directly  at 
the  fuel  cell  cathode  by  flowing  the  oxygen-rich  carrier  fluid  through  the  fuel  cell  electrode. 
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Inis  allows  the  fuel  cell  to  remain  at  ambient  pressure  while  eliminating  the  inefficient 
concentration  or  compression  of  the  extracted  oxygen.  g  incmcient 

T^tsdirect  earner  feed  fuel  cell  (CFFC)  is  a  new  development  in  DARPA's  artificial  gill 

program  which  re-introduces  the  possibility  of  using  fuel  cells  for  underwater  power 
generation.  r 

The  exact  chemical  reaction  mechanism  is  not  yet  known,  but  the  overall  reaction  is: 

(LM)2  -  02  +  2H20  +  4e"  --->  2LM  +  40H' 


where 


(LM)2  -  02  is  oxygenated  carrier,  and 
LM  is  deoxygenated  carrier 


To  counteract  the  OH*  generated,  electrochemical  oxidation  of  hydrogen  is  employed  on  the 
other  side  of  a  cation-exchange  membrane.  The  hydrogen  ion  generated  is  transported 
across  the  membrane  and  neutralizes  the  OH". 

The  reaction  at  the  anode  is, 

2H2  —  >  4H+  +  4e" 

The  hybrid  fuel  cell  design  is  a  new  concept  which  will  require  further  development  in 
order  to  increase  its  overall  efficiency.  Details  of  the  design  can  be  found  in  Appendix  C. 

ALWATT  HYDROGEN  GENERATOR 

Hydrogen,  of  course,  is  the  other  requirement  for  the  fuel  cell  to  produce  power.  Although 
there  arc  many  possible  sources  of  hydrogen,  the  next  section  will  show  that  the  Alwatt 
system  provides  a  much  safer  and  volumetrically  superior  option  for  the  CFFC.  The 
Alwatt  is  a  corrosion  cell  where  the  aluminum  is  reacted  with  seawater  to  provide  some 
electrical  power  and  a  significant  amount  of  hydrogen  as  a  by-product.  This  permits  the 
use  of  the  reactant  in-situ  to  provide  a  source  of  gaseous  hydrogen  at  ambient  pressure. 

The  overall  reaction  is: 

A1  +  3H20  — >  A1(0H)3  + 1 H2 

More  details  about  the  Alwatt  system  can  be  found  in  Appendix  D. 

The  following  are  some  salient  features  of  the  CFFC  system  designed  to  provide  a  power 
of  10  watts  continuously  for  a  duration  of  one  year. 


4 


Long  Endurance  Underwater  Power  System 
Aquanautics  Corporation _ 


Quarterly  Report 
April  1989 -June  1989 


Weight  (kg) 


Fuel  Cell  2.2  3.33 

Gill  2.5  2.2 

Alwatt  45.4  113.0 

Seawater  Pump  *  * 

Carrier  Pump  +  Control  *  * 


Total  50.1  118.53 


♦Included  in  Alwatt  volume  and  weight. 


Power  Generated 
/Consumed  (Watts) 

+10.0 

0 

+2.30 

-1.15 

-1.15 


+10.00 


2.2  Design 

2.2.1  System  Energy  Balance 

The  chemical  reaction  taking  place  in  the  Alwatt  is  as  follows: 

2A1  +  6  H20  — ~>  2  Al(OH)3  +  3  H2 

From  this  equation,  we  can  derive  the  following: 

2  moles  of  Aluminum  give  3  moles  of  Hydrogen, 

Or,  1  g.  of  Aluminum  gives  1.2444  liters  of  Hydrogen. 

From  the  Alupower  report1  on  the  performance  of  the  Alwatt,  we  know  that 

1  g.  of  Aluminum  gives  1.5  A-hrs.  at  0.25  V, 

1  g.  of  Aluminum  gives  0.25  X  1.5  =  0.375  W-hrs  of  energy. 

Therefore,  the  energy  density  of  the  Alwatt  is  0.375  W-hrs/g  or  225  W-hrs/kg.  of 
aluminum 

To  calculate  the  energy  derived  from  the  Hydrogen  evolved  by  the  Alwatt, 

1  mole  of  Hydrogen  gives  2  X  96,500  C  or  A-sec,  (2  electrons/molecule  of  H2), 

I  mole  of  Hydrogen  gives  193,000  A-sec  or  53.61 1 1  A-hrs, 

I I  of  Hydrogen  gives  2.3934  A-hrs,  (since  1  mole  =  22.4  liters),  and 

1  g.  of  Aluminum  gives  2.9784  A-hrs  (since  1  g  of  Aluminum  =  1.2444 1  of  H2). 

If  the  Hydrogen  can  be  oxidized  in  a  fuel  cell  at  a  potential  of  0.6  V  (which  is  half  the 
theoretical  voltage  that  can  be  obtained  from  that  process),  we  can  derive  the  following: 

1  g.  of  Aluminum  gives  2.9784  A-hrs  X  0.6  V  =  1.787  W-hrs  of  energy  via  H2, 
The  energy  density  of  the  Hydrogen  fuel  cell  is  1.787  W-hrs/g  or  1787  W-hrs/kg  of 
Aluminum. 


1B.M.L.  Rao,  W.  Halliop,  Alupower  Report  on  Hydrogen  and  Contaminant  Study,  submitted  to 
Aquanautics,  Reported  to  DARPA,  Quarterly  Report  of  Oct  '89  -  Dec.  '89. 
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When  the  two  systems  are  superimposed,  the  energy  density  is  2162  W-hrs/kg  of 
Aluminum,  which  is  roughly  a  factor  of  5.75  higher  than  the  energy  density  of  the 
Aluminum  alone  in  the  Alwatt,  where  the  Hydrogen  by-product  is  wasted. 

The  above  calculations  are  based  on  a  unit  quantity  of  aluminum  and  do  not  relate  to  system 
volume  or  weight.  The  assumptions  were: 


1)  For  Alwatt: 


Operating  Voltage  =  0,25  V 

A-hr./g  of  A1  a  1.5  (theoretical  3) 

These  above  numbers  were  taken  from  Alupower  report,  Alupower  being 
manufacturer  of  the  Alwatt  system. 

2)  For  Aquanautics  CFFC  System: 


Fuel  Cell  Voltage  «  0.6 

Hydrogen  Utilization  «  60% 

The  fuel  cell  voltage  constitutes  a  performance  target,  which  will  be  discussed  later. 
But  utilizing  the  full  amount  of  hydrogen  generated  in  Alwatt  is  virtually 
impossible,  The  major  loss  will  be: 


1 )  Through  the  fuel  cell  membrane  by  diffusion 

2)  In  the  seawater  flowing  through  Alwatt  by  dissolution 

The  magnitudes  of  these  losses  are  calculated  in  Section  3.1. 

From  these  calculations  we  can  conclude  that  the  amount  of  H2  needed/watt 
assuming  a  fuel  cell  voltage  of  0.6V)  is: 


1)  For  fuel  cell  reaction  ■  11.6ml/min 

2)  For  loss  through  fuel  cell  =  5  ml/min 

3)  For  loss  through  seawater  through  Alwatt  »  1.7  ml/min 


TOTAL  -  18.3  ml/min 

Aluminum  needed  ■  14.7  mg/min 

From  the  Alupower  report,  it  was  determined  that  the  Alwatt  at  6000  m  will  have  a  starting 
voltage  of  around  0.22  volts.  However,  at  0.22  volts  and  0.15  A-hr/g  of  Al,  14.7  mg/min 
will  generate  around  290  mW.  Since  the  voltage  will  drop  as  Al  gets  consumed  and  due  to 
shunt  current  losses,  a  further  20%  drop  was  assumed,  leading  to  230  mW, 

2.2.2  Aquanautics  System  Comparison  With  Lithium  Thionyl  Chloride 
Battery 

In  the  preceding  section,  various  types  of  oxygen  sources,  hydrogen  sources,  and  energy 
sources  have  been  evaluated.  In  each  case,  the  technology  offering  the  smallest  size  has 
been  selected  and  the  Aquanautics  CFFC  system  defined  as  an  integration  of  these  best 
options.  This  section  deals  with  the  evaluation  and  sizing  of  the  integrated  system  and  its 
comparison  with  systems  based  on  competing  technologies. 
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The  size  of  the  Aquanautics  CFFC  system  was  calculated  by  first  computing  individual 
component  size.  TTie  equations  and  theory  behind  the  sizing  of  individual  components  are 
discussed  in  the  succeeding  subsections. 

Table  2.1  shows  a  summary  spreadsheet  indicating  all  the  important  parameter  values  of  the 
system  for  a  10W,  1  year  system. 

2.2.3  Calculation  of  the  Gill  Volume  and  Weight 

The  gill  volume  and  weight  were  calculated  on  the  basis  of  the  following  two  equations  that 
were  derived  to  describe  the  seawater  pumping  power  and  the  oxygen  flux  through  the  gill: 


P 

Q3  \2  p 


13-5093^&^)0-2'ln^V'2 


=  4.3354  [ 


D2  \r2  i 


$<$-l)2  Q2  do5 


(1/3) 


(2.2.1) 

(2.2.2) 


where,  P  is  the  theoretical  power  needed  for  seawater  pumping  (W), 

V  is  the  volume  of  the  gill  fiber  bed  (m^), 

Q  is  the  seawater  flow  rate  through  the  gill  (m3/s), 

1  is  the  thickness  of  the  gill  fiber  bed  (m), 

p  is  the  density  of  the  seawater  (kg/m3), 

4  is  the  fiber  spacing  factor,  such  that  x  do  is  the  distance  between 
the  centers  of  adjacent  fibers  in  one  layer, 
u  is  the  kinematic  viscosity  of  seawater  (m^/s), 

D  is  the  diffusivity  of  oxygen  through  seawater  (m^/s), 

Csi  is  the  initial  concentration  of  oxygen  in  seawater  (std  m3/m3), 
N  is  the  oxygen  flux  desired  for  the  system  (std  m^/s),  and 
do  is  the  outer  diameter  of  the  gill  fibers  (m). 


The  derivation  of  these  two  governing  equations  is  shown  in  Section  4. 

The  net  power  output  of  the  system  was  used  to  specify  the  value  for  the  seawater  pumping 
power  P  in  Equation  2.2.1.  For  example  if  the  net  power  was  assumed  to  be  10  W,  for  a 
system  of  that  size,  it  was  determined  that  the  Alwatt  would  produce  2.3  W.  Based  on  the 
assumption  that  the  Alwatt  power  would  be  used  for  the  various  hotel  loads,  1.15  W  or 
50%  of  the  Alwatt  output  was  set  aside  for  seawater  pumping.  A  seawater  pumping 
efficiency  of  10%  was  assumed  and  the  theoretical  value  of  the  pumping  power  P  was 
assumed  to  be  0.115  W.  The  oxygen  flux  desired  was  calculated  from  the  net  output 
power,  which  was  assumed  to  be  supplied  by  the  Hydrogen-Carrier  fuel  cell.  The  net 
output  voltage  of  the  fuel  cell  was  assumed  and  the  cell  current  calculated  from  the 
equation, 


(2.2.3) 
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where,  Pn  is  the  net  power  output  (W),  and 
E  is  the  cell  potential  (V). 

The  oxygen  flux  needed  was  calculated  from  the  current  value  obtained  from  Equation 
2.2.3  by  using  the  following  equation: 

N  =  jpj  (2.2.4) 

where,  n  is  the  electron  count  (4  electrons/molecule  of  oxygen), 

F  is  the  Faraday  Constant  (96,500  C/mole),  and 
N  is  the  oxygen  flux  (mol/s). 
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Table  2.1 


19-Jul-89 

Model  of  Direct  Feed  Fuel  Cell  and  Gill 

With  Carrier  Directly  Reduced  at  the  Fuel  Cell  Cathode 


IS  yin 

"Onus - 

1000  W 

W 

Net  Power 

1.08  Eo 

V 

Fuel  CeUEO 

0.48  Ep 

V 

Losses  in  Fuel  Cell  Efficiency 

300  i 

A/raA2 

Fuel  Cell  Current  Density 

25  Rf 

1/a 

Fuel  Cell  Packing  Density 

1  L 

year 

Endurance 

1500  rho  .c 

kg/m*3 

Density  of  the  fuel  cell 

2500  rho_a 

kg/mA3 

Density  of  Alwan  Cell 

Outputs 

Units 

5.56E+00  Af 

mA2 

Fuel  Cell  Area 

4.96E-02  Qs 

mA3/s 

Seawater  Flow  Rate 

0.60  E 

volts 

Fuel  Cell  Voltage 

2.22E-01  Vf 

mA3 

Fuel  Cell  Volume 

1.67E-01  Vg 

mA3 

Theoretical  GUI  Volume 

50%  ng 

Gill  volume  overhead  allowance 

2.50E  01  Vga 

aA3 

Actual  Gill  Volume 

4.72E-01  Vfg 

mA3 

Gill  Volume  +  Fuel  Cell  Volume 

1.82E+00  Vh 

mA3 

Theoretical  Hydrogen  Source  (Alwatt)  Volume 

40%  nh 

Packing  efficiency  +  H2  generation  efficiency  of  Alwan 

4.54E+00  Vha 

mA3 

Actual  volume  of  Hydrogen  source  (Alwatt) 

5.01E+00  Vt 

01*3 

Total  Volume 

2.21E+02  Wg 

kg 

Gill  Weight 

1.13E+04  Wa 

kg 

Alwatt  Weight 

3.33E+02  Wc 

kg 

Fuel  CeU  Weight 

1.19E+04  Wt 

kg 

Total  Weight 

Calculations 

Units 

1666.67  If 

imp 

Fuel  Cell  Current 

9.67E-05  Qo 

mA3/s 

Oxygen  Consumption  rate 

1.93E  04  Qh 

aA3/s 

Hydrogen  consumption  rate 

8.40E-02  P/w  Jen 

W/kg 

Gravimetric  Power  density 

2.00E+02  P/v_den 

W/mA3 

Volumetric  Power  density 

7.36E+02  E/w.den 

W-hrAcg 

Gravimetric  Energy  density 

1.75E+06  E/vJen 

W-hrymA3 

Volumetric  Energy  density 
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Gill  V  &  Q*  f(P) 


Inputs 


- F 

- rTTW - — 

wTTTS37W 

L 

0.025  m 

x,b  3019.45049 

rho 

1025.8  k*/mA3 

y.c  0.46829456 

z 

4 

A  0.03223949 

u 

1.762  IE-06  mA2/s 

B  7.78125408 

D 

2.1E-09  mA2/s 

Cs 

0.003  mA3/mA3 

N 

9.6718E-05  mA3/s 

do 

0.0003  m 

Outputs 

Q  0.04958748  mA3/s 


Q  f  of  Q  fp  of  q  Qp  expBQA(2/3)  Q2thirds 

0.04958748  7.7139E-13  2.55857433  0.04958748  2.85839971  0.13497335 
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Conversion  of  Gill  Volume  to  Weight 
Input  Parameters 


di 

0.00025  m 

do 

0.0003  m 

V 

0.16654814  mA3 

den_o 

2200  kg/mA3 

den_f 

1200  kg/m*3 

den_c 

1140  kg/mA3 

zeta 

4 

Calculations 

Vo 

0.08327407  mA3 

Vf 

0.00999219  mA3 

Vc 

0.02270953  mA3 

Output  Gill  Weight 

W 

221.082461  kg. 
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A  spreadsheet,  "Gill  V  &  Q  =  f(P),"  was  used  to  compute  the  gill  volume  V  and  seawater 
flow  rate  Q  from  Equations  2.2.1  and  2.2.2  using  the  Newton-Raphson  method.  Once  the 
gill  volume  V  was  calculated  another  spreadsheet,  "Gill  Volume  and  Weight,"  was  used  to 
compute  the  weight  of  the  gill.  The  gill  volume  V  obtained  from  Equations  2.2.1  and  2.2.2 
represent  only  the  volume  of  the  fiber  bed.  This  value  of  V  was  multiplied  by  a  factor  of 
1.5  to  obtain  a  total  gill  volume  including  the  support  structures  and  the  seawater  pump. 
Therefore,  the  gill  packaging  efficiency  was  assumed  to  be  67%.  The  overall  gill  weight 
was  now  assumed  to  consist  of  three  essential  components  as  shown  in  Equation  2.3.5. 

Wg  =  W0  +  Wf  +  Wc  (2.2.5) 

where,  Wg  is  the  total  weight  of  the  gill  (kg), 

W0  is  the  weight  of  the  gill  support  structure  and  pump  (kg), 

Wf  is  the  weight  of  the  gill  fibers  (kg),  and 
Wc  is  the  weight  of  the  carrier  in  the  gill  (kg). 


The  weights  of  the  individual  gill  components  were  obtained  by  simply  multiplying  the 
volumes  associated  with  those  components  by  a  density.  The  gill  support  structure  volume 
was,  as  mentioned  before,  0.5  V.  The  fiber  volume  and  carrier  volume  were  computed  to 


be. 


Vf  = 


PJV 

4$ 


U 


d‘2> 

d?} 


(2.2.6) 


Vc  = 


p  V  dj2 
4^  do2 


(2.2.7) 


where,  Vf  is  the  volume  occupied  by  the  fibers  (m3), 

Vc  is  the  volume  occupied  by  the  carrier  in  the  fibers  (m3),  and 
dj  is  the  inner  diameter  of  the  fibers  (m). 


The  densities  that  were  assumed  for  the  various  components  are  as  follows: 


po  =  2200  kg/m3  (gill  packaging  and  pump), 

pf  =  1200  kg/m3  (gill  fiber  material,  polypropylene  or  polyethylene),  and 
pc  =  1 140  kg/m3  (carrier). 

The  weight  of  the  gill  was  therefore  calculated  from  the  equation, 


Wg  =  Vcpc  +  Vfpf  +  VqPq  (2.2.8) 

2.2.4  Calculation  of  Fuel  Cell  Volume  and  Weight 

The  fuel  cell  volume  and  weight  were  calculated  on  the  basis  of  the  expected  performance 
characteristics.  The  current  calculated  from  Equation  2.2.3  was  used  to  size  the  fuel  cell. 
The  current  density  was  assumed  to  be  30  mA/cm2  or  300  A/m2.  Current  densities  two 
orders  of  magnitude  higher  than  this  have  been  achieved  in  conventional  fuel  cells  using 
Hydrogen  and  Oxygen  and  therefore,  30  mA/cm2  represents  a  rather  conservative  estimate. 
This  current  density  was  divided  into  the  total  cell  current  to  yield  the  cell  area  as  follows: 
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Ac  = 


I 

i 


(2.2.9) 


where,  Ac  is  the  area  of  the  fuel  cell  membrane  (m2), 

I  is  the  fuel  cell  current  (A),  and 
i  is  the  current  density  (A/m^). 

The  fuel  cell  packaging  efficiency  was  calculated  to  be  25  m^/m-  of  volume.  This  indicates 
that  one  cubic  meter  of  fuel  cell  volume  was  assumed  to  hold  25  m^  of  membrane  area. 
This  number  corresponds  to  the  packing  efficiency  of  the  multipurpose  ABC-250  cell  that 
has  been  developed  at  Aquanautics  for  its  oxygen  extraction  system.  The  same  cell  can  be 
employed  as  a  fuel  cell  by  incorporating  the  appropriate  electrode  configuration  and  by 
using  the  fuel  cell  membrane. 

Therefore  the  volume  of  the  fuel  cell  Vfc  was  calculated  using  Equation  2.2.2. 

Vfc  =  ^  (2.2.10) 

From  the  ABC-250  experience,  it  is  known  that  the  cell  density,  which  is  a  composite  of 
that  of  the  cell  parts,  is  about  1500  kg/m^.  The  same  number  was  used  for  the  purposes  of 

sizing  the  fuel  cell  also,  given  the  adaptability  of  the  ABC-250  design.  This  number  is 

conservative  given  the  fact  that  a  rcjgular  electrochemical  cell  has  liquid  electrolyte  in  both 
chambers  (anode  and  cathode)  while  the  fuel  cell  application  would  entail  use  of  liquid 
electrolyte  -  carrier  in  this  case  -  only  in  the  cathode  chamber.  Multiplying  the  fuel  cell 
volume  Vfc  by  the  cell  density  yielded  the  cell  weight  Wfc  as  shown  in  Equation  2.2.1 1, 

Wfc  =  1500  Vfc  (2.2.11) 

2.2.5  Calculation  of  the  Alwatt  Volume  and  Weight 

Of  the  three  major  components  of  the  Aquanautics  system  for  oxygen  extraction  and  power 
generation,  the  Alwatt  is  the  only  one  whose  size  depends  on  the  endurance  of  the  system. 
The  Alwatt  system  was  sized  based  on  the  quantity  of  Hydrogen  needed  to  produce  the 
required  W-hrs  of  energy.  The  capacity  of  the  system  was  calculated  from  Equation 
2.2.12. 


C  = 


Pn  T 
E 


(2.2.12) 


where,  C  is  the  capacity  of  the  system  (A-hrs), 

T  is  the  system  endurance  (hrs), 

Pn  is  the  net  power  output  of  the  system  (W),  and 
E  is  the  cell  potential  of  an  individual  fuel  cell  (V). 

The  amount  of  Hydrogen  required  over  the  lifetime  of  the  system  was  calculated  from  the 
current  capacity  of  the  system  as  follows: 
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Qh 


3600  C 
n  F 


(2.2.13) 


where,  Qh  is  the  quantity  of  Hydrogen  storage  required  (moles) 

n  is  the  electron  count  for  Hydrogen  (2  electrons/molecule  H2),  and 
F  is  the  Faraday  Constant  (96,500  C/mole). 

The  equation  for  die  reaction  taking  place  in  the  Alwatt  is  as  follows: 

2A1  +  6H2O  - — >  2  Al(OH)3  +  3H2- 

From  this  chemical  reaction,  it  is  obvious  that  to  produce  1  mole  of  Hydrogen,  2/3  moles 
or  18  g.  of  Aluminum  are  required.  Therefore  the  quantity  of  Aluminum  required  to  store 
Qh  moles  of  Hydrogen  is  calculated  by  multiplying  Qh  into  18.  The  Alwatt  volume  could 
then  be  calculated  by  using  a  density  for  the  Aluminum  of  2800  kg/m^.  In  addition,  a  few 
other  factors  had  to  be  taken  into  account  to  compute  the  Alwatt  volume.  They  were: 

1)  Loss  of  Hydrogen  through  diffusion  and  dissolution  in  seawater, 

2;  Packaging  inefficiencies  in  the  Alwatt  cell,  and 
3)  Support  structures  for  the  Alwatt 

Accordingly  the  theoretical  volume  for  the  Alwatt  derived  in  the  preceding  paragraph  was 
multiplied  by  a  factor  of  2.5  to  account  for  the  above  contingencies  and  the  final  Alwati 
volume  obtained.  A  density  of  2500  kg/m 3  was  attributed  to  the  Alwatt  and  the  weight 
calculated  by  multiplying  this  density  into  the  volume  of  the  Alwatt.  The  density  of  the 
Alwatt  is  lower  than  the  density  of  the  Aluminum  itself,  because  of  voids  between  plates  of 
Aluminum  fuel  and  the  fact  that  the  packaging  material  into  which  the  Alwatt  is  built  is 
likely  to  be  non-metallic  and  made  of  plastic  material. 

The  volume  of  the  carrier  pump  and  the  controls  were  not  included  in  this  analysis,  and 
they  are  expected  to  be  negligible  when  compared  to  the  size  of  the  three  major  components 
analyzed  here.  Moreover,  some  of  the  packaging  efficiencies  assumed  in  this  study  are 
generous  enough  that  they  could  be  assumed  to  include  allowances  for  the  carrier  pump  and 
control  electronics,  if  any.  For  example,  the  volume  of  the  carrier  pump,  for  such  a 
system,  from  laboratory  experience  is  not  likely  to  be  more  than  0.75  liters.  The  carrier 
pump  will  be  pumping  anywhere  between  200  and  600  ml/min  of  carrier  fluid  through  the 
system  at  current  performance  levels  and  pumps  to  do  this  are  available  commercially  in 
small  packages.  It  is  therefore  assumed  that  employing  conservative  packaging  efficiencies 
for  the  three  major  components  would  provide  enough  contingency  for  the  carrier  pump 
and  system  controls.  This  is  not  expected  to  lead  to  significant  errors. 

2.2.6  Calculation  of  Competing  System  Sizes 

The  sizes  of  competing  systems  were  calculated  based  on  published  energy  densities.  The 
specifications  of  most  power  systems  include  data  on  the  gravimetric  and  volumetric  energy 
density  and  this  data  was  used  in  the  following  equations  to  determine  the  system  weight 
and  volume: 


Ws 


PnT 

cSg 


(2.2.14) 
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Vs  =  (2.2.15) 

e$v 

where,  Ws  is  the  weight  of  the  competing  system  (kg). 

Vs  is  the  volume  of  the  competing  system  (m^), 

Pn  is  the  power  output  of  the  system  (W), 

T  is  the  endurance  of  the  system  (hr) 

e$g  is  the  gravimetric  energy  density  of  the  system  (W  hr/kg),  and 
eSv  is  the  volumetric  energy  density  of  the  system  (W  hr/m^). 

Implementation  of  System  Sizing  Model 

The  sizing  model  discussed  in  the  above  subsections  was  implemented  in  spreadsheets. 
The  gill  sizing  was  accomplished  by  solving  Equations  2.2.1  and  2.2.2  iteratively  by  using 
the  Newton-Raphson  method.  To  do  this,  the  two  equations  were  combined  to  yield  an 
equation  of  the  form, 

f(Q)  =  BQ@/3)  +  In  (1  -  ^)  =  0  (2.2.16) 

where,  A  and  B  are  constants  which  combine  the  parameters  D,  u,  r,  C$i, 

N,  P,  do,  and  1, 

Q  is  the  seawater  flow  rate. 

Equation  2.2.16  was  then  solved  using  the  Newton-Raphson  method  and  the  gill  volume 
determined  by  substituting  the  value  of  Q,  thus  determined,  into  Equation  2.2.2.  The 
iterative  process  was  implemented  in  a  spreadsheet  called  "Gill  V  &  Q  =  f(P)"  and  the 
output  namely  the  values  of  V  and  Q  were  used  by  the  other  linked  spreadsheets.  The  gill 
volume  was  converted  into  a  gill  weight  by  the  spreadsheet  "Gill  Volume  and  Weight”, 
while  the  volume  and  weight  of  the  fuel  cell  and  the  Alwatt  were  computed  in  the  main 
spreadsheet  called  "Deep  Power  Model".  The  main  spreadsheet  also  served  as  the  point  at 
which  the  assumptions  to  the  model  were  introduced  and  the  results  of  the  model  supplied 
to  tables  for  tabulation  of  the  system  size  as  a  function  of  endurance  or  fuel  cell  voltage. 

The  equations  describing  the  size  of  competing  systems  were  implemented  in  the  main 
spreadsheet  "Deep  Power  Model"  also.  This  was  done  so  that  the  tables  used  for  plotting 
would  have  the  sizes  of  both  the  Aquanautics  CFFC  as  well  as  the  alternatives.  Printouts 
of  the  spreadsheets  are  attached  for  reference. 
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3.0  FUEL  CELL 


3.1  Calculation  of  Hydrogen  Leak  Through  Fuel  Cell  and  Seawater  Flow 
Through  Alwatt  System 

3.1.0  Background 

The  DARPA  project  calls  for  a  study  into  the  feasibility  of  extracting  oxygen  from  seawater 
and  combining  this  oxygen  with  hydrogen  in  a  fuel  cell  to  produce  electrical  power.  The 
hydrogen  is  generated  from  an  aluminum- seawater  cell  as  a  by-product.  The  system  is 
being  considered  for  application  at  depths  of  up  to  6000  m  under  the  ocean  and  is  being 
designed  to  operate  at  die  ambient  pressure.  This  implies  that  the  hydrogen  in  the  fuel  cell 
would  be  at  pressures  of  up  to  600  atmospheres.  The  high  hydrogen  pressure  in  the  fuel 
cell  creates  a  problem,  namely  one  of  diffusion  to  the  sea  through  the  fuel  cell  membrane, 
carrier  and  gill  membrane  to  the  seawater.  Normally,  since  the  atmosphere  does  not 
contain  significant  quantities  of  hydrogen,  the  hydrogen  content  of  seawater  can  be 
considered  to  be  zero.  This  causes  existence  of  a  steep  hydrogen  concentration  difference 
between  the  fuel  cell  anode  and  the  seawater  surrounding  the  system.  Since  the  carrier  and 
the  seawater  come  in  contact  with  one  another  through  the  gill,  there  is  a  potential  for 
hydrogen  diffusion  through  the  gill  to  the  ocean.  This  diffusive  flux  has  to  be  estimated  so 
that  appropriate  measures  can  be  enacted  to  minimize  the  loss  of  hydrogen  to  the  seawater. 

3.2.0  Objective 

The  objective  of  this  study  is  to  determine  the  rate  of  hydrogen  diffusion  from  the  fuel  cell 
to  the  seawater.  This  is  to  be  accomplished  by  considering  the  hydrogen  permeability  of 
the  fuel  cell  membrane  as  well  as  the  mass  transport  coefficients  for  the  carrier  and  seawater 
streams.  Another  objective  of  this  exercise  is  to  set  up  a  model  for  the  diffusion  of 
hydrogen  through  the  system  so  that  remedial  measures  can  be  identified.  Also  such  a 
model  is  expected  to  help  identify  the  system  component  that  controls  the  flux  of  hydrogen 
through  the  system. 

3.3.0  Model  for  Hydrogen  Transport 


3.3.1  Overall  Transport  Equation 


Figure  3.3.1  Diffusion  Path  for  Hydrogen  from  Fuel  Cell  to  Seawater 
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The  system  schematic  shown  in  Figure  3.3.1  identifies  the  various  components  of  the  mass 
transfer  coefficient  for  hydrogen  flux  from  the  fuel  cell  to  seawater. 

The  hydrogen  transport  problem  can  be  modelled  using  an  electrical  analog  with  resistors  in 
series.  The  effective  mass  transport  coefficient  for  hydrogen  flux  from  the  fuel  cell  to  the 
seawater  can  therefore  be  defined  by  Equation  3.3.1  which  is  normalized  with  respect  to 
the  area  of  the  fuel  cell  membrane. 


1  _  1  1  Ac  Ac  Ac  Ac  Aq 

ke  knP  kcP  Agjkgc  Qc  +  Ag  k$  +  Qs  +  Ag  kg 


(3.3.1) 


where,  ke  is  the  effective  mass  transport  coefficient  for  hydrogen  transport 
from  the  fuel  cell  to  the  seawater  (m/s), 
kn  is  the  hydrogen  permeability  of  the  membrane  (m/s), 
p  is  the  fraction  of  the  total  membrane  area  exposed  by  the  fuel  cell 
anode, 

Ac  is  the  area  of  the  membrane  in  the  fuel  cell  (m2), 

kc  is  the  mass  transport  coefficient  for  the  earner  in  the  cathode  of 

the  fuel  cell  (m/s), 

kgC  is  the  mass  transport  coefficient  for  the  carrier  stream  in  the  gill 
(m/s), 

Agj  is  the  surface  area  presented  to  the  carrier  by  the  inner  surfaces 
of  the  gill  hollow  fibers  (m2), 

ks  is  the  mass  transport  coefficient  for  the  seawater  stream  flowing 
through  the  gill  (m/s), 

Ag  is  the  area  presented  to  the  seawater  by  the  outside  surface  of  the 
gill  hollow  fibers  (m2), 

kg  is  the  hydrogen  permeability  of  the  gill  membrane  (m/s), 

Qc  is  the  flow  rate  of  the  carrier  through  the  system  (m^/s),  and 
Qs  is  the  flow  rate  of  seawater  through  the  gill  (m^/s). 


The  permeabilities  of  the  fuel  cell  membrane,  kc,  and  the  gill  membrane,  kg,  are  usually 
specified  in  terms  of  a  flux  rate  per  unit  area  per  unit  partial  pressure  difference 
(mol/s/m2/atm).  This  can  be  converted  to  the  standard  mass  transfer  coefficient  units  of 
m/s  by  using  Henry's  Law.  The  solubility  of  hydrogen  (given  in  mol/m^/atm)  can  be 
divided  into  the  published  permeability  to  yield  the  familiar  units  of  m/s.  What  is 
happening  here  is  that  hydrogen  is  diffusing  from  a  gaseous  phase  to  a  dissolved  phase 
where  the  solvent  is  a  liquid.  This  means  that  the  concentrations  are  usually  expressed  in 
two  different  units,  atm  for  the  gas  phase  and  mol/m^  for  the  dissolved  state.  Using 
Henry’s  Law  allows  conversion  of  the  membrane  permeability  referred  in  terms  of  gas 
phase  concentrations  (atm)  to  liquid  solution  phase  concentrations  (mol/m^). 

Another  point  to  be  noted  with  respect  to  Equation  3.3.1  is  that  the  overall  mass  transport 
coefficient  is  being  referred  in  terms  of  the  fuel  cell  membrane  area.  This  is  because  the 
fuel  cell  membrane  was  determined  to  be  the  crucial  and  controlling  component  of  the 
hydrogen  diffusion  problem.  Therefore,  optimization  of  the  system  will  proceed  by 
adjusting  the  fuel  cell  area  and  analyzing  its  implications. 
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The  magnitude  of  the  hydrogen  transport  can  be  gaged  from  the  value  of  the  term  ke,  and 
the  relative  importance  of  its  various  constituents  determined  from  the  values  of  the 
individual  components  in  Equation  3.3.1.  The  hydrogen  transport  can  now  be 
characterized  by  the  equation. 

Qh  =  keAcAC,  (3.3.2) 

where,  Qh  is  the  flux  rate  of  hydrogen  (mol/s  or  m3  @  STP/s),  and 

AC  is  the  difference  in  hydrogen  concentration  between  the  fuel  cell 
anode  and  the  seawater  (mo  1/m 3  or  @  STP/m^). 

For  the  sake  of  clarity,  moles  are  used  as  the  preferred  units  forexpressing  hydrogen  flux. 
Since  the  expression  for  ke  had  been  referred  to  the  fuel  cell  area  in  Equation  3.3.1,  use  of 
Ac  in  Equation  3.3.2  is  permitted.  It  can  therefore  be  seen  from  Equation  3.3.2  that  in 
order  to  reduce  the  magnitude  of  Qh,  the  value  of  Ac  has  to  be  reduced.  Since  the  units  for 
ke  are  m/s,  the  units  for  AC  should  be  mol/m 3.  The  concentration  of  hydrogen  in  the  gas 
phase  at  the  fuel  cell  anode  (given  in  atmospheres  of  pressure)  will,  therefore  have  to  be 
converted  to  solution  concentrations  (mol/rrP)  by  using  Henry’s  Law.  The  fuel  cell  anode 
hydrogen  partial  pressure  is  thus  converted  into  a  solution  concentration  that  would  exist  in 
equilibrium  with  that  partial  pressure. 

3.3.2  Determination  of  Fuel  Cell  Area 

The  fuel  cell  area  is  the  most  crucial  component  of  the  hydrogen  diffusion  problem.  The 
area  of  the  fuel  cell  membrane  is  determined  from  the  current  density  that  is  found  to  be 
ideal  for  operation.  The  current  density  at  the  operating  point  is  determined  from  the  output 
power  requirements  and  the  cell  voltage  at  that  current  density.  The  cell  voltage  is  generally 
of  the  form. 


Vc  =  V0-ai, 


(3.3.3) 


where,  Vc  is  the  fuel  cell  voltage  under  load  (V), 

V0  is  the  fuel  cell  voltage  under  no  load  conditions  (V), 
i  is  the  current  density  (A/m^),  and 
a  is  a  constant  (Cl  -  m2). 

The  value  of  the  constant  a  is  dependent  on  a  number  of  factors,  such  as  mass  transfer, 
activation  kinetics  and  the  cell  impedance.  One  component  of  the  cell  impedance  is  the 
resistance  of  the  membrane.  This  resistance  is  linearly  dependent  on  the  thickness  of  the 
membrane  and  therefore  this  can  be  easily  incorporates  into  Equation  3.3.3  to  give. 


Vc=  V0  (p  +  trm)i, 


(3.3.4) 


where,  P  is  a  constant  (Q  -  m2), 

t  is  the  membrane  thickness  (m),  and 

Tm  is  the  specific  resistance  of  the  membrane  (Cl  - 
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This  equation  is  derived  by  splitting  the  constant  a  from  Equation  3.3.3  into  two 
components,  one  related  to  the  membrane  resistance  and  the  other  encompassing  the 
remaining  components  of  the  fuel  cell  losses.  This  helps  to  identify  another  parameter  that 
afreets  Hydrogen  diffusion,  namely  the  thickness  of  the  fuel  cell  membrane. 

If  the  desired  output  power  can  be  denoted  by  P,  the  cell  area  can  be  calculated  from  the 
equation  relating  cell  voltage  and  current  to  power, 

Ac  =  (3.3.5) 


where,  P  is  the  gross  power  output  of  the  fuel  cell  (W),  and 

Ac  is  the  fuel  cell  area  (m^). 


At  this  stage,  it  is  necessary  to  divide  the  total  cell  area  into  a  number  of  individual  cells  to 
enable  determination  of  the  mass  transport  coefficient  on  the  carrier  side  of  the  cell  or  in  its 
cathode  compartment.  Usually  the  cell  dimension^  are  restricted  by  designs  that  are 
commercially  available,  and  it  can  be  assumed  here  that  the  available  cells  have  a  length  of  1 
and  a  breadth  of  b.  Therefore  area  of  each  individual  cell  is  given  by, 

ac  =  bl  (3.3.6) 

where,  ac  is  the  area  of  the  individual  cell  (m^), 

b  is  the  breadth  of  the  cell  (m),  and 
1  is  the  length  of  the  cell  in  the  direction  of  carrier  flow  (m). 

The  number  of  cells  is  therefore  given  by 


nc  = 


Ac 

«c 


(3.3.7) 


where,  nc  is  the  number  of  cells 


3. 1.3. 3  Membrane  Permeability 

The  membrane  permeability  can  now  be  calculated  given  the  membrane  thickness  t. 
Manufacturers  of  membranes  that  are  commonly  used  in  fuel  cells  publish  permeabilities 
for  membranes  of  unit  thickness.  An  example  would  be  the  units  used  to  specify  the 
hydrogen  permeability  of  Nafion  which  is  given  in  mol  -  mil/s/m^/atm .  The  actual 
permeability  of  the  membrane  can  be  determined  from  this  value,  as  well  as  its  thickness 
and  the  solubility  of  Hydrogen  in  the  liquid  medium  to  which  the  transport  is  occurring. 
This  yields, 


*n  * 


Kn 
Sh  t 


(3.3.8) 


where,  Kn  is  the  specific  permeability  of  the  cell  membrane  (mol-mil  /s 
An^/atm),  and 

Sh  is  the  solubility  of  hydrogen  in  seawater  (raol/m^/atm). 
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3.3.4  Carrier  Flow  Rate 


The  carrier  flow  rate  determines  two  constituents  of  the  total  hydrogen  flux,  namely 
diffusion  on  the  cathode  side  of  the  fuel  cell  membrane  and  the  capacity  for  transmission 
through  the  canier  itself.  To  determine  the  carrier  flow  rate,  it  is  necessary  to  consider  the 
hydrogen  and  oxygen  consumption  rates.  The  hydrogen  consumption  rate  is  directly 
dependent  on  the  fuel  cell  current  and  is  given  by  the  equation 


Qhc  = 


i  Ac 

TTf 


(3.:.9) 


where,  Qhc  is  the  rate  at  which  hydrogen  is  consumed  in  the  fuel  cell 

(mol/s), 

n  is  the  number  of  electrons  involved  in  the  oxidation  of  each 
hydrogen  molecule  (2).  and 
F  is  the  Faraday  constant  (96,500  C/mol). 


The  hydrogen  consumption  rate  yields  the  oxygen  consumption  rate  by  stoichiometry. 


Qoc  =  0.5  Qhc. 


(3.3.10) 


where,  Qoc  is  the  rate  at  which  oxygen  is  consumed  (mol/s). 

The  carrier  flow  rate  is  also  dependent  on  the  fractional  conversion  of  carrier  in  the  fuel  cell 
and  the  concentration  of  carrier  in  solution.  The  fractional  conversion  is  a  term  that  defines 
the  proportion  of  earner  that  is  dcoxygenated  in  the  fuel  cell.  Hence, 

<3-31» 


where,  Qc  is  the  canier  flow  rate  (m^/s). 

Me  is  the  molar  concentration  of  the  carrier  solution  (mol/m^),  and 
fc  is  the  fractional  conversion  of  the  carrier  in  the  fuel  cell. 

The  use  of  the  multiplier  2  in  Equation  3.3.1 1  is  indicative  of  the  fact  that  the  carrier  is 
dimeric,  namely  that  it  takes  two  molecules  of  the  carrier  to  bind  one  oxygen  molecule.  To 
calculate  the  mass  transfer  coefficient  on  the  carrier  side  of  the  cell,  the  flow  rate  of  carrier 
through  each  cell  is  calculated  as  per  the  equation, 

qc  =  ^  (3.3.12) 


where,  qc  is  the  flow  rate  of  carrier  per  cell  (m^/s). 

3.3.5  Mass  Transport  of  Hydrogen  to  the  Carrier 

The  mass  transport  coefficient  for  the  transport  of  hydrogen  to  the  carrier  in  the  cathode 
compartment  of  the  fuel  cell  can  be  calculated  based  on  the  Leveque  correlation  for  mass 
transport  into  flow  in  conduits.  This  correlation  is  expressed  as: 
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Sh  =  1.62  Re0  33  ScO*33 

where,  Sh  is  the  Sherwood  Number, 

Re  is  the  Reynolds  Number, 

Sc  is  the  Schmidt  Number, 

df  is  the  hydraulic  diameter  of  the  cathode  conduit  (m),  and 
1  is  the  length  of  the  flow  path  in  the  cathode  compartment  (m). 

The  equation  describing  the  transport  of  hydrogen  to  the  carrier  in  the  fuel  cell  cathode 
compartment  does  not  account  for  the  likely  presence  of  a  high  surface  area  electrode  in  that 
compartment.  The  presence  of  such  an  electrode,  which  may  take  on  the  form  of  a  carbon 
felt,  will  serve  to  increase  the  resistance  to  hydrogen  flux.  The  equation  listed  above, 
therefore,  represents  a  worst  case  scenario.  The  dimensionless  numbers  referred  to  above 
are  defined  below. 

Re=  icjcdf 
lie  b  tf 

where,  tf  is  the  thickness  of  the  cathode  compartment  (m), 

rc  is  the  density  of  the  carrier  (kg/m3),  and 
me  is  the  dynamic  viscosity  of  the  carrier  (N  s/m^). 


Sc 


Ik 

rc  Dc 


where. 


Dc  is  the  diffusivity  of  hydrogen  in  carrier  (m^/s). 


Sh  = 


kc  df 

Ik 


where,  kc  is  the  mass  transport  coefficient  for  the  carrier  side  of  the  cell 
(m/s). 


The  hydraulic  diameter  of  the  conduit  df  is  defined  in  traditional  terms  as  four  times  the 
cross-sectional  area  divided  by  the  perimeter,  and  is  therefore  represented  by. 


df 


4  b  tf 


Using  the  above  equations,  the  carrier  side  mass  transport  coefficient  was  determined  to  be 


(3.3.13) 
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3.3.6  Gill  Sizing 

Gill  sizing  is  important  to  the  hydrogen  diffusion  problem  because  hydrogen  exits  the 
system  through  the  gill.  It  is,  however,  not  likely  that  the  gill  will  control  the  hydrogen 
flux  because  of  its  large  interfacial  area.  It  is  likely  not  to  provide  significant  resistance  to 
hydrogen  transport.  The  gill  resistance  to  hydrogen  transport  is  composed  of  three 
components  just  as  in  the  case  of  oxygen  transport  into  the  system.  These  three 
components  are  the  carrier  side  boundary  layer,  the  membrane  permeability  and  the 
seawater  side  boundary  layer.  To  calculate  the  effects  of  the  gill  operation  on  the  hydrogen 
flux,  it  is  necessary  to  size  the  gill  first.  The  gill  sizing  is  dependent  on  the  oxygen  flux 
that  is  desired  for  system  power  generation.  The  gill  design  equations  yield,  the  volume  of 
the  gill  Vg,  the  area  of  the  gill  Ag,  and  the  length  of  the  fibers  a.  Also  as  a  result  of  gill 
sizing,  the  flow  rate  of  seawater  can  be  calculated  given  the  extraction  efficiency  of  the 
system.  The  gill  design  also  yields  a  thickness  for  the  gill  which  is  denoted  here  by  lg. 
The  length  of  the  fibers  in  a  gill  of  square  cross-section  can  be  determined  by  the  equation. 


(3.3.14) 


where,  Vg,  lg  and  a  are  defined  above. 

The  number  of  fibers,  which  is  important  in  determining  the  flow  and  transport 
characteristics  through  each  fiber,  can  be  calculated  from  the  surface  area  of  the  gill  and  the 
outer  diameter  of  each  fiber. 


"f  *  (3-3.15) 

r  d0  a 

where,  nf  is  the  number  of  gill  fibers,  and 

do  is  the  outer  diameter  of  the  gill  fibers  (m). 

The  flow  rate  of  carrier  through  each  gill  fiber  can  therefore  be  calculated  using  the 
equation, 

Qcf  =  ~  (3.3.16) 

where,  qcf  is  the  carrier  flow  rate  through  each  fiber  (m3/s). 

3.3.7  Mass  Transport  Coefficients  in  the  Gill 

The  mass  transport  characteristics  for  hydrogen  diffusion  through  the  gill  are  determined 
using  the  correlations  of  non-dimensional  numbers  that  have  been  derived  for  various  flow 
situations  and  geometries.  For  the  flow  of  carrier  through  the  fibers,  the  Leveque 
correlation  shown  in  the  beginning  of  Section  3.3.5  is  employed.  This  gives  an  expression 
for  the  mass  transport  coefficient  shown  below 


(3.3.17) 
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where,  kgC  is  die  mass  transport  coefficient  for  the  earner  side  in  the  gill 

(m/s), 

dj  is  the  inner  diameter  of  the  gill  hollow  fiber  (m),  and 

qcf,  Dc.  Pc.  Pc.  and  a  have  all  been  defined  in  Sections  3.3.5  and 
E.3.7. 

For  the  hydrogen  flux  characteristics  on  the  seawater  side,  a  correlation  derived  by  Yang 
and  Cussler  is  used.  This  correlation  was  derived  for  flow  perpendicular  to  a  hollow  fiber 
bed.  This  correlation  relates  the  Reynolds,  Schmidt  and  Sherwood  numbers  in  the  familiar 
fashion, 

Sh  =  1.38  Re°  34  Sc®  33 

Using  the  usual  definitions  of  the  above  dimensionless  numbers,  the  seawater  side  mass 
transfer  coefficient  if  derived  to  be 


where,  us  is  the  seawater  velocity  through  the  gill  (m/s) , 

ps  is  the  seawater  density  (kg/m3), 

Ps  is  the  dynamic  viscosity  of  the  seawater  (N  s/m?-), 

Ds  is  the  diffusivity  of  hydrogen  jirough  seawater  (m2/s), 
ks  is  the  mass  transfer  coefficient  for  hydrogen  diffusion  on  the 
seawater  side  (m/s),  and 

do  is  the  outer  diameter  of  the  gill  hollow  fibers  (m). 

The  seawater  velocity  term  us  has  to  be  determined  from  the  seawater  flow  rate,  taking  into 
account  the  fact  that  the  entire  gill  cross-section  (a  x  a)  is  not  available  for  seawater  flow. 
This  term  was  derived  for  a  gill  consisting  of  fibers  (£  do)  apart  and  arranged  in  layers 
where  adjacent  layers  are  perpendicular  to  each  other.  The  seawater  velocity  is  derived  to 
be, 


us 


QA1 

(£-l)2  a2 


(3.3.19) 


where,  £  is  the  spacing  factor  for  the  fibers  in  the  plane  perpendicular  to  the 

seawater  flow  direction,  such  that  the  center  to  center  distance 
between  adjacent  fibers  is  ^  do- 

The  mass  transfer  coefficients  thus  determined  can  be  entered  into  Equation  3.3.1  to 
determine  the  effective  mass  transfer  coefficient  and  the  hydrogen  flux  from  the  fuel  cell  to 
the  seawater. 
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3.4.0  Application  of  the  Model  to  a  Case  Study 

The  model  developed  based  on  the  equations  derived  in  Section  3.3.0  was  implemented  in 
a  spreadsheet.  This  was  done  to  examine  the  Hydrogen  diffusion  characteristics  and 
evaluate  the  effect  of  various  fuel  ceil  operating  parameters  on  the  extent  of  the  Hydrogen 
loss  through  the  fuel  cell  membrane.  The  parametric  analysis  was  performed  by  plotting 
the  fraction  of  the  total  Hydrogen  consumed,  that  was  lost  to  diffusion  at  various  depths. 

3.4.1  Solubility  of  Hydrogen  in  Seawater 

Initially,  to  examine  the  Hydrogen  diffusion  phenomenon,  data  was  collected  for  the 
solubility  of  Hydrogen  in  seawater,  since  it  is  clear  that  the  Hydrogen  diffusion  seeks  to 
establish  such  an  equilibrium.  As  part  of  the  feasibility  study  undertaken  by  Alupower2  the 
solubility  of  Hydrogen  in  seawater  was  calculated  based  on  published  data  for  Hydrogen 
solubility  in  pure  water.  The  results  of  this  study  are  summarized  in  Table  3.4.1.  Included 
in  Table  3.4.1  are  the  solubility  of  Hydrogen  in  seawater  as  a  function  of  depth  and  a 
conversion  of  this  solubility  into  a  Henry's  Law  Constant.  The  Henry's  Law  Constant  is 
usually  defined  in  the  equation, 

P  *  (3.4.1) 

where.  p  is  the  partial  pressure  of  the  gas  under  consideration  (atm), 

%  is  the  mole  fraction  of  the  gas  dissolved  in  the  solvent  (moles 
solute/mole  solution),  and 

H  is  the  Henry's  Law  Constant  (atm  mole  solution/mole  solute). 


It  is  obvious  here  that  the  Henry's  I^w  Constant  depicted  in  Table  3.4.1  is  the  reciprocal  of 
the  constant  used  in  Equation  3.4.1.  The  constant  in  Table  3.4.1  has  units  that  allow  it  to 
be  used  as  the  parameter  Sh  in  Equation  3.3.8. 


TABLE  3.4.1 

DEPTH 

PRESSURE 

SOLUBILITY 

HENRY'S  CONST 

(m) 

(atm) 

(std  m^/nv*) 

(std  m^/rn^/atm) 

0 

1.0000 

0.019 

0.0190 

500 

50.6575 

0.93 

0.0184 

1000 

100.3151 

1.855 

0.0185 

3000 

298.9452 

5.347 

0.0179 

5000 

497.5753 

8.569 

0.0172 

7000 

696.2054 

11.646 

0.0167 

10000 

994.1506 

15.679 

0.0158 

3.4.2  Assumptions  in  the  Case  Study 


This  section  lists  the  main  assumptions  made  in  evaluating  the  Hydrogen  diffusion 
problem.  These  assumptions  arc  divided  into  appropriate  sections  dealing  with  individual 
components  as  below: 

Reference  already  cited. 
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1)  Power  System  Assumptions: 

System  Power  =  1000  W 

Parasitic  Power  Consumption  (Pumps  &  Controls)  =  230  W 
Seawater  Pump  Power  Consumption  =  1 15  W 

2)  Fuel  Cell  Assumptions:  (Data  from  Nafion  Specification  Sheet  with  Giner 

Electrode) 


Fuel  Cell  Membrane  Permeability  (Kn)  = 

Fuel  Cell  Membrane  Porosity  (p)  - 

Fuel  CeU  Sizing  (ABC-250)  (1 X  b) 

Fuel  Cell  Membrane  Specific  Resistance  (rm)  = 

Slope  of  Cell  Potential  -  Current  Density  curve  (P)  = 

Thickness  of  Fuel  Cell  Chamber  (tf)  = 


520  cm3-mil/(ft2  hr  atm) 
30% 

lOcmX  25  cm 
0.1463  Qm2/m 

0.001148  Qm2/rn 
7.5  mm 


3)  Carrier  Property  Assumptions: 

Diffusivity  of  Hydrogen  through  carrier  (Dc) 
Density  of  carrier  (pc) 

Dynamic  Viscosity  of  carrier  (pc) 
Concentration  of  Carrier  in  Solution  (Mc) 


=  5.36  X  10-9  m2/s 
=  1 140  kg/m  3 
*  1.455  cP 

=  0.4  M  or  400  moles/m^ 


4)  Gill  Sizing  Assumptions: 

Gill  Thickness  (lg) 

Fiber  Outer  Diameter  (do) 

Fiber  Inner  Diameter  (di) 

Fiber  Spacing  Factor  ($) 

5)  Seawater  Property  Assumptions: 

Initial  Oxygen  Concentration  in  Seawater  (Cs0 
Density  of  Seawater  (ps) 

Kinematic  Viscosity  of  Seawater  (us) 
Diffusivity  of  Oxygen  through  Seawater  (D0s) 
Diffusivity  of  Hydrogen  through  Seawater  (D§) 

3.4.3  Parametric  Variations  in  Case  Study: 


0.05  m 
400  pm 
350  pm 

4 


3  ml  @  STP/1 
1025.8  kg/m3 
1.7621  X  10~6  m2/ 
2.1  X  10-9  m2/s 
5.36  X  IQ'9  m2/s 


The  parametric  analysis  was  conducted,  as  mentioned  earlier,  to  evaluate  the  fraction  of 
Hydrogen  consumption  that  is  lost  through  diffusion  to  seawater  at  vanuur  depths.  This 
analysis  was  conducted  by  varying  the  values  of  four  parameters  that  were  thought  to 
impact  the  diffusion  rate  directly.  These  parameters  are: 
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1)  Fuel  Cell  Membrane  Thickness  (t), 

2)  Fuel  Cell  Current  Density  (i), 

3)  Fractional  Conversion  in  the  Fuel  Cell  (fc),  and 

4)  Open  Circuit  Cell  Fuel  Cell  Potential  (VG). 

Table  3.4.2  shows  the  values  of  the  various  parameters  that  were  used  in  the  parametric 
analysis.  This  table  also  depicts  the  baseline  values  used  in  the  analysis.  The  baseline 
value  represents  the  value  of  the  parameter  that  was  held  constant  while  the  variation  of 
another  parameter  was  being  studied.  For  instance,  in  studying  the  effect  of  the  membrane 
thickness  by  varying  its  value,  the  value  of  the  current  density  was  held  constant  at  the 
baseline  value,  namely  30  mA/cm^. 


Table  3.4.2 

Memb.  Thick. 

Current  Density 

Fractional  Conv 

OC  Cell  Poten. 

Comments 

(mils) 

(mA/cm^) 

(%) 

(V) 

5.0 

10 

2.5 

0.75 

7.5 

20 

10.0 

0.85 

10.0 

30 

5.0 

0.95 

Baseline 

12.5 

40 

15.0 

1.05 

15.0 

50 

25.0 

1.15 

3.4.4  Method  of  Application  of  Model 

The  model  was  applied  to  evaluate  the  different  values  of  the  mass  transport  coefficients  in 
Equation  3.3.1.  For  the  purposes  of  this  analysis,  it  was  assumed  that  the  gill  membrane 
permeability  kg  did  not  control  the  Hydrogen  mass  transport  and  was  therefore  ignored. 
The  other  mass  transport  coefficients  calculated  included  the  following: 

kn  Membrane  Permeability 

kc  Fuel  Cell  Cathode  Mass  Transfer  Coefficient 

kgC  Carrier  Side  Mass  Transfer  Coefficient  for  the  gill 

ks  Seawater  Side  Mass  Transfer  Coefficient  for  the  gill 

Qc/Ac  Carrier  Capacity 

Qs/Ac  Seawater  Capacity 

To  perform  these  calculations,  three  intcr-linked  spreadsheets  were  devised.  The  first 
spreadsheet  contained  the  information  regarding  the  solubility  of  Hydrogen  at  various 
depths  that  is  shown  in  Table  3.4.1.  This  table  was  used  as  a  look-up  street,  and  was 
referred  to  whenever  information  regarding  Hydrogen  dissolution  equilibrium  was  needed. 
The  second  sheet  was  devoted  to  gill  sizing.  Since  the  Hydrogen  diffusion  is  affected  by 
gill  sizing  and  the  seawater  flow  rates,  this  information  was  considered  necessary.  The  gill 
sizing  spreadsheet  was  set  up  to  receive  the  seawater  pumping  power  information  from  the 
main  sheet,  which  constituted  the  actual  calculation  of  the  hydrogen  diffusion  mass 
transport  coefficients. 

The  power  output  of  the  system  was  entered  through  the  main  spreadsheet  and  an  allotment 
for  seawater  pumping  power  calculated  from  the  net  power.  This  pumping  power 
information  was  shared  by  the  gill  sizing  spreadsheet,  which  used  it  along  with  tne  pill 
sizing  and  seawater  property  assumptions  to  calculated  a  gill  size.  Tne  gill  sizing 
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spreadsheet  used  the  two  simultaneous  equations  derived  in  Section  4.  These  equations 
were  solved  for  the  given  power  value  to  yield  gill  volume  and  the  seawater  flow  rate.  The 
gill  sizing  spreadsheet  used  the  Newton-Raphson  method  to  solve  the  gill  design  equations. 

The  gill  sizing  information  and  the  seawater  flow  rate  from  the  gill  sizing  spreadsheet  were 
then  used  by  die  main  spreadsheet  to  calculate  the  various  mass  transport  coefficients  as  per 
the  equations  developed  in  Section  3.3.0. 

3.5.0  Results  and  Discussion  of  the  Parametric  Analysis 

As  mentioned  above,  the  effects  of  the  various  parameters  are  investigated  by  plotting  the 
percentage  of  Hydrogen  lost  to  diffusion  as  a  function  of  depth.  The  effect  of  each 
parameter  is  investigated  in  a  separate  section. 

3.5.1  Effect  of  Fuel  Cell  Membrane  Thickness 


Figure  3.5. 1  Hydrogen  Diffusion  for  various  Fuel  Cell  Membrane  Thickness  values 

Figure  3.5. 1  depicts  the  variation  of  Hydrogen  loss  due  to  diffusion  as  a  function  of  depth 
for  various  values  of  the  fuel  cell  membrane.  The  values  in  the  legend  of  Figure  3.5.1  have 
units  of  mils  (milli-inchcs).  As  expected,  the  diffusion  loss  is  least  for  the  thickest 
membranes  and  the  increase  of  the  diffusion  rate  as  a  function  of  depth  is  not  linear.  This 
is  because  the  solubility  curve  for  Hydrogen  is  not  linear.  A  point  to  note  here  is  that 
increased  membrane  thickness  will  cause  an  increase  in  the  ohmic  losses  in  the  cell.  This  is 
due  to  the  increased  electrical  resistance  of  thicker  membranes.  This  drop  in  voltage  would 
normally  cause  the  need  for  higher  currents  to  maintain  the  net  power  output.  An  increased 
current  at  constant  current  densities  implies  greater  cell  areas,  which  would  in  normal 
circumstances  cause  increased  hydrogen  diffusion.  In  the  model  however,  the  contribution 
of  the  term  trra  -  which  represents  the  resistance  of  the  membrane  -  is  smaller  than  that  of  b 
which  represents  other  contributions  to  losses  in  the  cell.  The  output  voltage  is  hence  not 
verv  highly  sensitive  to  the  membrane  thickness.  Therefore  the  negative  effect  of  Increased 
cell  losses  are  more  than  offset  by  the  beneficial  impact  of  decrease  in  fuel  cell  membrane 
permeability. 
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3.5.2  Effect  of  Fuel  Cell  Current  Density 

Figure  3.5.2  depicts  the  variation  of  the  diffusion  loss  of  Hydrogen  for  various  values  of 
the  current  density  (units  in  the  legend  are  A/m^)  in  the  fuel  cell.  As  expected  the  current 
density  has  a  strong  effect  on  the  Effusion  loss.  This  is  because  the  area  of  the  fuel  cell  is 
inversely  proportional  to  the  current  density.  Lower  current  densities  therefore  necessitate 
the  use  of  larger  fuel  cell  areas  causing  the  potential  for  greater  diffusion  of  Hydrogen.  The 
increased  diffusion  rates  at  low  current  densities  are  slightly  offset  by  the  fact  that  fuel  cell 
ohmic  and  kinetic  losses  are  reduced  at  low  current  densities.  Therefore  the  cell  output 
voltage  is  increased  and  the  fuel  cell  current  needed  to  provide  the  necessary  power 
decreases.  This  beneficial  impact  is,  however,  not  very  significant  compared  to  the  impact 
of  the  increased  cell  area  as  is  evident  from  Figure  3.5.2. 


Figure  3.5.2  Hydrogen  Diffusion  for  various  Fuel  Cell  Current  Density  values 

3.5.3  Effect  of  Fractional  Conversion  of  Carrier  in  Fuel  Cell 

The  fractional  conversion  is  defined  as  that  fraction  of  the  total  number  of  moles  of  carrier 
presented  to  the  fuel  cell  that  is  converted.  Conversion,  here,  refers  to  the  electrochemical 
release  of  oxygen  for  purposes  of  producing  power  from  the  fuel  cell.  The  concept  of 
fractional  conversion  is  normally  tied  into  the  current  density,  that  is  higher  the  current 
density,  greater  the  number  of  carrier  molecules  being  reduced,  and  hence  greater  the 
fractional  conversion,  For  purposes  of  this  analysis,  the  fractional  conversion  and  the 
current  density  aurc  treated  separately  to  evaluate,  numerically,  the  effect  of  each.  In  a 
working  device,  it  is  expected  that  the  larger  fractional  conversions  would  be  obtained  at 
high  current  densities. 
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Figure  3.5.3  Hydrogen  Diffusion  for  various  Fractional  Conversion  values 

Figure  3.5.3  depicts  the  variation  of  the  Hydrogen  diffusion  with  different  values  of  the 
fractional  conversion.  As  expected  the  Hydrogen  loss  decreases  with  increasing  fractional 
conversion.  This  effect  of  higher  fractional  conversion  is  directly  manifested  in  smaller 
carrier  flow  rates.  This  causes  a  decrease  in  the  contribution  of  the  carrier  flow  to  the 
overall  mass  transport  coefficient  and  therefore  lower  Hydrogen  diffusion  rates.  When,  as 
is  expected  in  a  real  process,  the  increased  fractional  conversion  is  coupled  with  the  higher 
current  density,  the  decrease  in  Hydrogen  diffusion  will  be  more  dramatic.  This  is  because 
the  beneficial  effects  of  decreased  carrier  flow  rate  will  be  coupled  with  that  of  the  lower 
fuel  cell  area. 

3.6.0  Hydrogen  Loss  as  Dissolved  in  Seawater  through  Alwatt  System 

Alwatt  needs  certain  amount  of  seawater  flow  to  prevent  clogging  of  their  system.  This  is 
quoted  to  be  0.2  A-hr./ml  of  electrolyte  or  seawater  (Reference  cited). 

Calculation  of  hydrogen  loss  via  seawater  flew  through  Alwatt: 

1  g  of  A1  produces  1 .2444  liters  of  H2 

1  g.  of  A1  provides  1.5  A-hr. 

ml  of  seawater/g  of  A1  to  prevent  clogging: 

83  AW fficlwrolytc  '  7'5 11,1  °f  of  A1 

Therefore,  1.2444  liters  of  H2  (a  1  g  of  Al)  would  need  7.5  ml  of  seawater. 

At  6000  m,  H2  solubility  is  10  ml/ml  of  seawater  (Reference  cited). 
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Therefore,  if  we  provide  the  n  inimum  seawater  flow  of  7.5  ml  of  seawater/g  of  A1 
consumed,  we  would  lose  75  ml  of  H2  out  of  1 .24  liters  produced.  With  a  safety  factor  of 
65%,  the  hydrogen  loss  will  be  10%.  So  design  point  seawater  flow  is  12.4  ml/g  of  A1 
consumed. 

The  loss  is  very  strongly  dependent  on  the  pressure  or  depth  of  operation.  At  6000  m,  the 
maximum  loss  is  about  30%  of  total  amount  of  hydrogen  coming  to  fuel  cell  (see  Figures 
3.5.1  -  3.5.3,  this  appendix).  Therefore,  for  each  watt  generated  by  fuel  cell,  we  need: 


1) 

for  electrochemical  reaction: 

(  1  1  of  H2  =  2.39  A-hr) 

1 1.6  ml/min 

2) 

for  loss  through  fuel  cell 

1 1.6  X  y 

5  ml/min 

3) 

for  loss  through  seawater  thru  Alwatt 

(11.6  +  5)  x  0.1 

Total  H2  need 

_ VL  ml/min 

18.26  ml/min 

The  development  of  a  Hydrogen  loss  model  was  found  to  be  important  to  the  analysis  of 
the  whole  deep  submergence  power  system.  The  system  analysis  was  performed  to 
compare  the  energy  density  and  performance  of  the  Aquanautics  system  with  those  based 
on  competing  technologies.  In  order  to  size  the  Aquanautics  system,  it  is  necessary  to 
determine  the  size  of  the  Hydrogen  source,  the  Alwatt  battery  in  this  case.  To  determine 
the  size  of  the  Alwatt,  it  is  imperative  the  total  required  capacity  of  Hydrogen  be  computed, 
taking  into  account  not  only  the  amount  that  is  likely  to  be  consumed,  but  also  that  which  is 
likely  to  be  wasted  in  diffusive  losses.  This  model  helped  in  identifying  the  magnitude  of 
the  latter,  and  the  former  could  be  had  from  stoichiometry,  leading  in  turn  to  an  accurate 
and  reliable  estimate  for  the  size  of  the  Alwatt  as  well  as  the  integrated  Aquanautics  deep 
submergence  power  system. 
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4.0  GILL:  Gill  Sizing  vs.  Seawater  Pumping 
4.0.0  Introduction 

This  study  was  performed  as  part  of  a  project  evaluating  the  feasibility  of  using  oxygen 
dissolved  in  seawater  to  produce  power  under  the  sea.  The  power  source  being  evaluated 
would  use  the  extracted  oxygen  and  hydrogen  from  another  source  for  purposes  of 
producing  power  in  a  fuel  cell.  The  oxygen  extraction  is  performed  using  Aquanautics’ 
proprietary  oxygen  absorbing  carrier  solution.  This  carrier  solution  is  to  be  fed  into  the 
fuel  cell  directly  where  oxygen  is  chemically  reduced  in  the  cathode  accompanied  by 
oxidation  of  the  hydrogen  fuel  to  produce  power. 

Oxygen  extraction  from  the  seawater  necessitates  bringing  the  seawater  and  the  carrier  in 
close  proximity  with  one  another,  so  that  net  oxygen  transfer  can  take  place.  An  important 
characteristic  that  defines  the  choice  of  this  unit  is  the  need  to  keep  the  carrier  and  the 
seawater  physically  separate  since  mixing  would  dilute  the  aqueous  carrier  solution.  As  a 
consequence,  membrane  contactors  were  deemed  necessary.  Considerations  of  size  and 
weight  dictate  the  use  of  hollow  fiber  membranes  which  enable  packaging  of  large 
interfacial  areas  in  small  volumes.  The  structure  within  which  the  hollow  fiber  membranes 
are  implemented  is  called  a  gill  because  of  its  functional  similarity  to  the  breathing  apparatus 
of  fishes.  In  the  system  contemplated  here,  the  carrier  would  be  pumped  through  the 
lumen  or  the  insides  of  the  hollow  fibers,  while  the  seawater  would  flow  on  the  outside  of 
the  fibers. 


Figure  4.0, 1  Schematic  of  an  Undersea  Power  Source  using  Aquanautics'  Oxygen  Carrier 


The  system  schematic  for  such  a  power  system  is  depicted  in  Figure  4.0.1.  The  carrier 
solution  is  circulated  by  the  carrier  pump  between  the  gill  where  it  absorbs  oxygen  from  the 
seawater  and  the  fuel  cell  where  the  oxygen  is  used  in  the  generation  of  power.  The 
hydrogen  for  the  fuel  cell  is  evolved  from  an  Aluminum-seawater  battery  which  also 


Long  Endurance  Underwater  Power  System  Quarterly  Report 

Aquanautics  Corporation _ 1 _ April  1989  -  Junc.1989 


produces  power.  The  power  produced  by  the  fuel  cell  is  considered  to  be  the  net  power 
output  of  the  system.  The  power  produced  by  the  Aluminum-seawater  battery,  on  the  other 
hand,  is  used  for  the  hotel  loads  of  the  system,  which  can  be  broadly  divided  into  the 
following: 

1)  Carrier  pump, 

2)  Seawater  pump  in  the  gill,  and 

3)  Controls  and  Instrumentation. 

Two  considerations  enter  into  the  design  of  a  power  system  using  the  Aquanautics  oxygen 
extraction  technology.  The  first  consideration  has  to  do  with  weight  and  volume.  In  order 
to  be  competitive  with  other  systems,  it  is  important  that  the  Aquanautics  system  be 
smaller,  both  in  weight  and  volume,  when  compared  to  existing  technologies.  The  second 
consideration  is  that  the  power  consumed  by  the  various  hotel  loads  described  above  be 
minimized.  Because  of  the  relative  paucity  of  dissolved  oxygen  in  seawater  -  typical 
concentrations  range  from  0.5  to  6  standard  milliliters  per  liter  of  seawater  -  the  seawater 
pumping  power  is  expected  to  be  the  largest  hotel  load  component.  Therefore,  system 
optimization  has  to  proceed  by  examination  of  the  gill  volume  and  the  seawater  pumping 
power  and  investigation  of  the  relationship  between  the  two. 

The  purpose  of  this  study  is  to  describe  a  model  which  can  be  used  to  optimize  the  design 
of  the  gill  component  of  the  power  system.  This  model  quantifies  the  relationship  between 
the  seawater  pumping  power  and  the  gill  volume  required  to  deliver  a  given  oxygen  flux 
rate. 

4.1.0  Governing  Equations 

4.1.1  Description  of  Gill  and  Assumptions  for  the  Model 

It  has  been  mentioned  in  the  introduction  to  this  study,  that  volume  considerations  dictate 
the  use  of  hollow  fiber  contactors.  Hollow  fiber  contactors  are  employed  in  a  variety  of 
arrays  and  present  large  interfacial  areas  packaged  efficiently.  Mass  transfer  across  a 
hollow  fiber  membrane  is  dictated  by  the  equation, 

N  -  k  A  AC 

where,  N  is  the  flux  of  oxygen  (mol/s), 

k  is  the  overall  mass  transfer  coefficient  (m/s) 

A  is  the  intcrfacial  area  or  area  of  the  gill  membrane  (m2),  and 
AC  is  the  difference  in  oxygen  concentration  between  the  seawater 
and  the  carrier  (mol/m^), 

The  mass  transfer  coefficient  k  is  composed  of  three  components,  This  is  depicted  in 
Figure  4,1,1.  The  three  components  are: 

1 )  Carrier  side  mass  transfer  coefficient  (kc), 

2)  Membrane  permeability  (km),  and 

3)  Seawater  side  mass  transfer  coefficient  (ks). 

The  overall  mass  transfer  coefficient  can  be  expressed  as  a  function  of  its  components  by 
using  the  following  relationship: 
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Figure  4.1,1  Schematic  of  Oxygen  Diffusion  through  a  Hollow  Fiber 

Membrane 

Of  the  three  components  that  constitute  the  mass  transfer  coefficient  k,  the  membrane 
permeability  can  be  ignored  because  it  is  usually  larger  than  the  mass  transport  coefficient 
on  the  seawater  side,  This  is  because  diffusion  through  liquids  is  much  slower  than  it  is 
through  a  thin  membrane,  On  the  carrier  side,  the  concentration  of  dissolved  oxygen  is 
zero,  This  is  because  the  oxygen  that  diffuses  into  the  carrier  from  the  seawater  is 
immediately  consumed  in  a  chemical  reaction  where  the  carrier  binds  the  oxygen.  The 
system  operates  with  an  excess  of  deoxygenated  carrier  ensuring  that  the  oxygen  diffusing 
into  the  carrier  is  bound  at  all  times,  hence  maintaining  a  concentration  gradient  at  all  points 
in  the  gill.  The  carrier  side  mass  transport  coefficient  can  therefore  be  ignored,  since  it  is 
dictated  by  the  reaction  kinetics  of  carrier  oxygenation,  which  is  faster  than  the  diffusion  of 
oxygen  through  the  seawater.  Therefore,  for  purposes  of  engineering  the  artificial  gill,  it 
can  be  assumed  that  the  overall  mass  transport  coefficient  is  tne  same  as  the  seawater  side 
mass  transport  coefficient, 

4.1,2  Theoretical  Power-Volume  Relationship 

The  design  of  a  gill  uses  two  types  of  tools,  The  first,  derived  from  fluid  mechanics,  is  the 
concept  of  a  friction  factor,  the  second,  bask’  to  muss  transfer,  is  the  idea  of  u  mass 
transfer  coefficient,  The  equations  describing  the  two  phenomena  are  introduced  here  and 
synthesized  into  one  power-volume  relationship. 
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A  friction  factor  f  is  most  commonly  defined3  in  terms  of  a  drag  force  F: 


F  = 


(4.1.1) 


where,  f  is  a  friction  factor, 

A  is  a  characteristic  area  of  the  gill  (m^), 

p  is  the  density  of  the  seawater  (kg/m^),  and 

v  is  the  average  velocity  of  the  seawater  flow  through  the  gill  (nVs). 


What  this  characteristic  area  should  be  depends  on  the  gill  design.  In  the  gills  that  are  being 
considered  for  this  application,  sea  water  is  pumped  through  narrow  channels  between 
hollow  fibers.  In  this  case,  the  dominant  force  is  the  drag  in  these  channels,  and  the 
characteristic  area  is  conventionally  assumed  to  be  that  of  outer  surface  of  the  gill 
membrane  hollow  fibers. 


The  mass  transfer  coefficient  k,  the  second  key  tool,  is  defined  in  the  equation  for  mass 
transfer, 


N  =  k  A  AC  (4.1.2) 

where,  N  is  the  oxygen  flux  rate  (mol/s), 

k  is  the  mass  transfer  coefficient  for  the  gill  (m/s), 

A  is  the  surface  area  of  the  gill  (m^),  and 

AC  is  the  concentration  gradient  across  the  gill  (mol/rn^). 

In  this  case,  A  is  again  the  surface  area  of  the  gill  membranes.  The  concentration 
differences  AC,  which  is  often  called  a  "driving  force,"  is  in  this  case  the  actual  oxygen 
concentration  In  the  sea  water  minus  that  in  equilibrium  with  the  carrier  solution.  Since  the 
oxygen  diffusing  into  the  carrier  solution  is  bound  immediately  in  a  chemical  reaction,  the 
concentration  of  oxygen  in  the  carrier  can  be  assumed  to  be  zero,  and  the  value  of  the  term 
AC  becomes  equal  to  the  concentration  of  oxygen  in  the  seawater. 

These  two  ideas  can  be  combined  to  determine  the  power  required  to  pump  seawater 
sufficient  to  allow  for  a  given  amount  of  mass  transfer  N.  The  seawatei  pumping  power 
required  to  facilitate  unit  flux  of  oxygen  can  therefore  be  expressed  by  the  ratio 


where,  P  is  the  power  required  for  seawater  pumping  through  the  gill  (W), 

The  ratio  of  P  to  N  should  be  minimized  in  an  optimum  gill  design.  To  achieve  this, 
Equation  4.1.3  is  written  In  the  form  of  a  dimensionless  equation, 


3  Baumeiiier,  T,,  Avallanc,  R.A.,  and  BaumoUtor,  T„  111,  Marks'  Standard  Handbook  for  Mechanical 
Engineers,  McGraw-Hill,  New  York  (1978), 
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P  V2  >  ffvu  N  \ 
pi  2  Q3  J  "  ^kJ^QACj 


(4.1.4) 


where,  Q  is  the  volume  tr<  rate  of  flow  of  seawater  through  the  gill  (m^/s), 

V  is  the  volume  of  the  gill  (m^),  and 

1  is  the  thickness  of  the  gill  in  the  direction  of  the  seawater  flow. 


The  volumetric  flow  Q  is  given  by. 


Q  = 


(4.1.5) 


For  example,  in  the  case  of  a  shell  and  tube  exchanger  where  seawater  and  carrier  flow 
along  parallel  paths,  the  term  1  would  refer  to  the  length  of  the  fibers.  In  the  case  of  flow 
across  a  bed  of  fibers,  where  the  seawater  flows  in  a  direction  perpendicular  to  that  of  the 
carrier  flow,  the  term  1  would  indicate  the  thickness  of  the  fiber  bed. 

This  dimensionless  equation  contains  three  dimensionless  groups.  One  of  these  [fv/2k],  is 
the  ratio  of  the  friction  factor  to  the  Stanton  number.  This  is  also  called  the  ratio  of  j- 
factors.  For  gases  without  form  drag  and  in  turbulent  flow,  this  ratio  is  about  one,  an 
equality  called  the  Reynolds  analogy4.  For  liquids  without  form  drag  and  in  turbulent 
flow, 

ft  *  (e)223  <4'16) 

where,  u  is  the  kinematic  viscosity  of  the  seawater  (m2/s),  and 
D  is  the  diffusivity  of  oxygen  through  seawater  (m2/s). 

This  relation  is  called  the  Chilton-Coburn  analogy.  Similar  relations  are  used  in 
overwhelming  detail  for  the  design  of  compact  heat  exchangers.  These  results  are  largely 
restricted  to  turbulent  flow,  The  non-dimensional  group  [fv /2k]  has  to  be  calculated  for  the 
artificial  gills,  where  form  drag  and  laminar  flow  are  usually  expected. 

The  second  dimensionless  group  on  the  right  hand  side  of  Equation  4.1.4  is  the  oxygen 
extraction  efficiency  (N/QAC).  this  number  is  essentially  the  ratio  of  flux  of  oxygen  to  the 
rate  at  which  it  is  presented  to  the  gill.  Clearly,  when  the  driving  force,  which  is  given  by 
the  concentration  difference  AC,  goes  to  zero,  this  number  becomes  large,  and  the  power 
required  will  also  become  large. 

The  third  dimensionless  group  (PV2/pl2g3)  in  Equation  4,1.4  represents  the  power  to 
pump  a  certain  amount  of  seawater,  Q,  through  a  gill  of  volume  V  and  thickness  1.  This 
ratio  is  different  than  that  used  in  conventional  studies  of  mixing,  where  the  impeller  size 
and  rotation  are  more  directly  considered.  However,  its  physical  significance  is  similar;  it 
represents  the  ratio  of  power  to  flow  within  the  gill.  The  importance  of  this  and  the  other 
dimensionless  groups  will  be  clearer  from  the  application  of  the  general  equation  to  the 
specific  case  or  designing  a  gill,  This  process  is  discussed  in  the  succeeding  sections. 


4  Bennett,  C.O.,  and  Myers,  J,E„  Momentum  Heat  and  Mass  Transfer,  McGraw-Hill,  New  York  (1962), 
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4.1.3  Power-VoHume  Equation  for  Flow  Across  a  Hollow  Fiber  Bed 

Equation  4.1.4  from  the  previous  subsection  is  of  a  general  nature  and  the  evaluation  of  the 
non-dimensional  term  [fv/2k]  is  dependent  on  the  nature  of  the  gill  geometry  and  the  mass 
transport  characteristics  associated  with  the  geometry.  Here  the  gill  is  assumed  to  consist 
of  hollow  fibers  of  outer  diameter  do  that  are  laid  against  each  other,  with  adjacent  layers  of 
hollow  fibers  running  perpendicular  to  one  another  (see  Figure  4.1.2).  The  overall 
geometry  of  the  gill  is  such  that  the  seawater  flows  through  a  square  cross-section  of  area 
Af  and  bed  depth  1.  The  cross-section  of  the  gill  is  assumed  to  be  square  and  of  dimension 

a  by  a.  The  distance  between  adjacent  fibers  is  given  by  £do  where  £  is  a  non-dimensional 
factor  that  defines  the  gap  between  the  fibers.  The  volume  of  the  gill  is  assumed  to  be  V 
and  the  seawater  flow  rate  Q.  The  oxygen  extraction  rate  is  assumed  to  be  N. 

It  must  be  noted  here  that  the  velocity  term  v  in  Equation  4.1.3  as  well  as  in  the  general 
form  for  the  power  (Equation  4.1.4)  is  the  average  interstitial  seawater  velocity  through  the 
gill.  Because  the  entire  gill  cross  section  is  not  available  for  seawater  flow  due  to  the 
presence  of  the  hollow  fibers,  the  velocity  cannot  be  calculated  from  the  seawater  flow  rate 
Q  and  the  frontal  area  of  the  gill  Af  as  in  Equation  4.1.5,  Instead,  it  has  to  be  calculated 
from  the  actual  area  available  for  flow,  namely  Aa.  The  ratio  of  the  frontal  area  to  the 
available  area  is  dependent  on  the  fiber  spacing  factor  as  shown  below. 


6£  .  _ji_ 

A*  (VO2 

where,  £  is  the  fiber  spacing  factor. 


Figure  4. 1 .2  Preferred  geometry  for  Gill 
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The  velocity  is  therefore  given  by, 


where, 


QU2 
v  (4-02 


(4.1.7) 


Q  is  the  seawater  flow  rate  (m^/s), 

V  is  the  volume  of  the  gill  (m3), 

1  is  the  thickness  of  the  gill  and  the  length  of  the  seawater  flow  path 
(m),  and 

4  is  the  fiber  spacing  factor. 


When  the  expression  forv  given  in  Equation  4.1.5  is  replaced  with  the  expression  given  in 
Equation  4.1.7  the  general  expression  for  seawater  pumping  power  becomes, 


P  v2  _  r  i  rf  vi  r  N  -i 

Q3l2p  ‘  1($.1)4J  *-2  kJ  l-QACJ 


(4,1.8) 


The  friction  factor  f  is  a  function  of  the  Reynolds  number  and  this  relationship  has  been 
determined  empirically  for  flow  over  banks  of  tubes.  From  the  published  data  the  friction 
factor  for  flow  over  a  hollow  fiber  bed  such  as  the  one  shown  in  Figure  4.1.2  can  be 
defined  as5, 


50  _  50  v 
Re  v  d0 


(4.1.9) 


where,  u  is  the  kinematic  viscosity  of  the  seawater  (m^/s), 

f  is  the  friction  factor,  and 

v  is  the  velocity  of  the  seawater  through  the  gill  (m/s). 

For  flow  perpendicular  to  a  hollow  fiber,  Yang  and  Cussler^  have  derived  a  correlation  for 
the  mass  transfer  coefficient  k.  The  mass  transfer  coefficient  can  be  determined  from  the 
Sherwood  Number  Sh  which,  for  laminar  flow  (Re  <  2000)  of  a  liquid  (Sc  «■  1000) 
perpendicular  to  a  hollow  fiber  array,  is  given  by 


Sh  =  1.38  Re0  34  Sc0.33t  (4.1.10) 

where,  Re  is  the  Reynolds  Number  (v  do/o), 

Sc  is  the  Schmidt  Number  (WD),  and 
Sh  is  the  Sherwood  Number  (k  do/D). 


For  simplicity  the  exponent  on  the  Reynolds  Number  is  assumed  to  be  0.33  (introducing  an 
error  of  approximately  1-5%  in  this  case)  also  and  the  equation  for  k  derived  to  be, 


I  _ 

I  ®Perry,  R.H.,  and  Green,  D.W.,  Chemical  Engineers’  Handbook,  McGraw-Hill,  New  York  (1984). 

I  ®Yang,  M.C.,  and  Cussler,  E,L.,"Designing  Hollow  Fiber  Contactors,''  AlChE  J.,  32, 1 1  (1986) 
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k 


=  1.38 


v(1/3)D(2/3) 


(4.1.11) 


where,  v  is  the  seawater  velocity  as  described  in  Equation  4. 1 .7  (m/s), 
do  is  the  outer  diameter  of  the  gill  hollow  fiber  (m),  and 
D  is  the  diffusivity  of  oxygen  through  seawater  (ro^/s). 

The  expression  for  v  from  Equation  4.1.7  is  used  in  Equation  4.1.11  to  yield  a  new 
expression  for  the  mass  transfer  coefficient  shown  in  Equation  4.1. 12. 


S(2/3)  Q<  1/3)  1(1/3)  d(2/3> 

=  (Tow)  wmjm 


(4.1.12) 


The  outer  diameter  of  the  fiber  can  be  expressed  as  a  function  of  other  parameters  using 
Equation  4.1.2.  The  area  of  the  gill  A  is  expressed  in  terms  of  the  number  of  fibers  and  the 
diameter  of  fibers  using  the  expression 


A  =  it  do  2  nt 

where,  a  is  the  length  of  the  fibers  (m),  and 

nt  is  the  total  number  of  fibers  in  the  gill. 

The  total  number  of  fibers  in  the  gill  can  be  calculated  by  multiplying  the  total  number  of 
fiber  layers  (in  the  direction  of  seawater  flow)  by  the  number  of  fibers  per  layer.  Using  a 

center  to  center  distance  of  ^d0  for  fibers  in  any  given  layer  and  assuming  that  adjacent 
layers  lie  against  each  other,  the  total  number  of  fibers  is  found  to  be, 

a  1 

"l  =  $d02 

where,  1  is  the  thickness  of  the  gill  (m). 


Substituting  for  nt  in  the  expression  for  the  gill  surface  area  yields  the  equation 
A  _  ft  a2  1  _  ft  V 
^  do  4  do 

Using  this  expression  for  the  gill  surface  area  in  Equation  4.1.2  yields  an  expression  for 
the  fiber  outer  diameter  that  is  shown  in  Equation  4.1.13. 


.  n  k  V  AC 
do  =  — Z - 

t  N 


(4.1.13) 


where,  k  is  the  mass  transfer  coefficient  on  the  seawater  side  (m/s), 

AC  is  the  oxygen  concentration  difference  in  the  gill  (mol/m^),  and 
N  is  the  oxygen  flux  rate  (mol/s). 
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The  volume  of  die  gill,  it  must  be  recognized  here,  is  the  frontal  area  multiplied  by  the 
thickness  of  the  gill  and  hence, 


Using  the  expression  for  do  to  eliminate  it  from  Equation  4,1.12  yields  a  new  expression 
for  the  mass  transfer  coefficient  shown  in  Equation  4.1. 14, 


,  5<4/5>  0<I/5M(!/5)d«/S)N<2/5> 

k  -  0.7675  -y^  * 

Equation  4. 1.9  can  be  rewritten  as, 


(4,1.14) 


Eliminating  do  from  the  above  equation  (using  Equation  4.1.13)  and  rearranging  terms, 
yields  an  expression  for  the  non  dimensional  group  [fv  /  2k)  as  shown  in  Equation  4.1 . 15, 


7  958 

7  958  1?  V  AC 


(4.1.15) 


By  substituting  the  expression  for  k  from  Equation  4,1.14  into  the  right  hand  side  of 
Equr  >u  4,1.15,  a  usable  and  final  expression  for  [fv  /  2k]  is  derived  as  shown  in 
Equation  4,1.16. 


««<"  (p)  Pf 


(4.1.16) 


This  expression  can  now  be  substituted  into  the  general  equation  for  seawater  pumping 
power  (Equation  4.1.8)  tc  yield  Equation  4.1.17. 

JLVt-  .  13  J093 -ill- («)  f-U-W2 

Q3l2p  13'“93  (^1)3.2  W  (qII/  (QACJ 
4.1.4  Equation  for  Mass  Transfer  Across  a  Hollow  Fiber  Bed 


(4.1.17) 


The  power-volume  relationship  expressed  in  Equation  4.1,17  does  not  constitute  a  well- 
posed  problem  because  there  are  two  unknowns.  The  intention  of  this  modelling  process  is 
to  determine  the  power  needed  to  pump  seawater  through  the  gill  as  a  function  of  the  gill 

volume.  In  Equation  4.1.17,  the  values  of  the  parameters  D,  p,  and  o  arc  given.  These 
are  properties  of  the  seawater  and  are  dependent  primarily  on  the  salinity  and  the 
temperature  of  the  seawater.  The  value  of  the  oxygen  flux  N  depends  on  the  desired 

capacity  of  the  power  generating  unit.  Further,  the  value  of  1,  $,  and  do  can  be  assumed. 
These  represent  the  gill  sizing  parameters  thet  define  the  hollow  fiber  arrangement  and  form 
the  factors  around  which  a  gill  can  be  optimized  given  the  desired  oxygen  flux.  Thus 
Equation  4.1.17  has  three  unknown  quantities  namely  AC,  Q,  and  P.  The  concentration 
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difference  and  the  seawater  flow  rate  are,  however,  related  through  the  oxygen  flux 
equation  (Equation  4.1.2)  and  therefore  the  term  AC  can  be  expressed  in  terms  of  the 
seawater  flow  rate  Q.  This  leaves  the  two  unknowns,  namely  Q  and  P. 

In  order  to  obtain  the  power  required  for  seawater  pumping,  the  value  for  the  seawater  flow 
rate  sufficient  to  provide  the  desired  flux  has  to  be  calculated  independently.  Fortunately, 
this  is  possible  by  resolving  Equation  4.1,2  to  yield  a  second  equation  relating  Q  to  V,  The 
derivation  of  such  an  equation  would  yield  a  well  posed  problem  with  two  unknowns  and 
two  equations, 

The  concentration  difference  term  in  the  last  non-dimensionul  group  in  Equation  4, 1.17  is 
usually  of  a  log-mean  difference  form  given  by, 


AC  - 


where,  C$i  is  the  oxygen  concentration  in  the  seawater  at  the  gill  inlet, 

Cso  is  the  oxygen  concentration  in  the  seawater  at  the  gill  exit, 

Cci  is  the  oxygen  concentration  in  the  carrier  at  tile  gill  inlet,  and 
Cco  Is  the  oxygen  concentration  in  the  carrier  at  the  gill  exit. 

All  the  concentration  terms  in  Equation  4.1.18  are  expressed  in  terms  of  mols/m^,  The 
oxygen  concentration  in  the  carrier  can  be  assumed  to  be  zero  at  all  locations  in  the  gill, 
since  we  are  assuming  that  an  excess  of  carrier  is  being  supplied  to  the  gill.  This  means 
that  the  oxygen  that  enters  the  carrier  in  solution  through  the  membrane  is  immediately 
bound  chemically  by  the  deoxygenated  carrier  molecules  that  exist  in  abundance.  This 
assumption  permits  ignoring  the  carrier  side  mass  transport  coefficient  when  calculating  the 
seawater  pumping  power.  Therefore  the  terms  Cci  and  Cco  ire  eliminated  from  Equation 
4.1.18  to  give  a  new  expression  for  the  concentration  gradient, 


AC 


(4.1.19) 


The  terms  Csi  and  Cso  can  also  be  related  to  the  oxygen  flux  and  the  seawater  flow  rate  by 
the  equation, 


N  -  Q(Csi-Cso). 
which  when  rearranged  yields, 

Cso  “  Csi  *  {3.  (4.1.20) 

When  this  expression  for  Cso  is  substituted  into  Equation  4.1.19  and  the  term  AC  in 
Equation  4.1.17  replaced  with  the  expression  from  Equation  4,1.19,  the  new  power 
equation  becomes, 
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P  V2 


Q3  l2 


*  13.5093 


it 


(t-o3 


.2 


(b)  (q^ 


.2  /■,_  . _M.2 

m  jo 

-LjEillsI 


(4.1.21) 


A  new  version  of  the  mass  transport  equation  is  now  developed  to  relate  the  flux  N  with  the 
seawater  flow  rate  Q  for  the  given  gill  design  and  volume.  This  is  done  by  replacing  the 
terms  k,  AC,  and  A  in  Equation  4.1.2.  To  do  thitt  a  new  expression  for  the  surface  area  A 
is  determined  from  the  relations  for  A  and  nt  in  the  preceding  section. 


a  n  V 

A  «*  - - 

k  do 

An  expression  for  the  term  AC  is  determined  by  substituting  the  value  for  Cso  from 
Equation  4,1,20  in  Equation  4,1,19,  When  the  new  expressions  for  A  and  AC  are 
introduced  into  Equation  4, 1 ,2  along  with  the  expression  for  k  from  liquation  4. 1,12  a  new 
expression  for  the  flux  N  results  as  shown  in  Equation  4.1,22, 


N 


1.38  ™ 


„D(2/3)v(2/3)  i(l/3) 


vw^/55  “oWBdjsBj . 

c.iT 


(4.1.22) 


By  eliminating  the  quantity  N  front  both  sides  of  Equation  4.1.22  an  expression  for  the  lust 
non-dimensional  group  of  Equation  4.1.21  can  be  determined  as  shown  in  Equation 
4,1.23. 


m{ 


“‘Sri 

"IS. 


»  4.3354 


$U/3)  ($-i)(2/3)  ~Q(3737do(5/35 


(4.1.23) 


The  value  for  Q  is  independently  determined  from  Equation  4.1.23,  given  the  geometry  of 
the  gill  (represented  by  V,  do,  1  and  $),  in  situ  properties  of  the  seawater  (represented  by 
C«i  and  D)  and  the  flux  desired,  N,  This  value  for  the  seawater  flow  rate  Q  is  then 
introduced  into  Equation  4.1.21  and  the  value  for  P  calculated.  A  parametric  analysis  can 

be  performed  by  changing  the  input  values  of  the  gill  sizing  parameters,  V,  do,  1.  and  £  to 
identify  the  optimum  gill  configuration.  Equations  4  '  ,21  and  4,1.23  therefore  form  a  well 
posed  problem  where  the  two  unknowns  Q  and  P  are  evaluated  as  functions  of  V. 

4.2.0  Application  of  the  Gill  Power*Volume  Model 

The  gill  power-volume  model  was  applied  and  a  parametric  analysis  performed  for  a  one 
standard  liter  per  minute  gill.  The  study  was  intended  to  identity  the  effect  of  the 
parameters  that  define  a  gill  design  on  the  power  to  pump  the  seawater  through  the  gill. 
The  various  values  for  which  the  model  was  applied  are  shown  in  Table  4.2.1,  and  the 
baseline  values  for  the  various  parameters  have  been  indicated.  For  example,  the  variation 

of  P  as  a  function  of  V  was  studied  for  four  different  values  of  £  keeping  the  other 
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parameters  constant  at  their  baseline  values,  This  process  was  repeated  for  each  individual 
parameter,  1,  do.  and  Csi* 


TABLE  4.2.1 

$ 

C.i 

Comments 

(microns) 

msm 

— 

rnn— 

Baseline 

2.5 

4(X) 

3,0 

1.50 

7.5 

500 

4,0 

2,25 

10.0 

600 

5.0 

4.50 

The  effects  of  the  various  parameters  on  the  power  for  seawater  pumping  are  presented  in 
sections  below,  Along  with  the  power  for  seawater  pumping,  two  other  parameters  namely 
seawater  flow  rate  and  the  pressure  drop  through  the  gill  were  also  evaluated,  The 
seawater  flow  rate  required  for  a  given  oxygen  flux  was  computed  from  Equation  4.1.23 
for  various  values  of  gill  volume,  and  the  value  for  the  seawater  pumping  power  calculated 
by  substituting  this  value  for  0  into  Equation  4.1,21.  The  value  for  the  pressure  drop 
through  the  gill  was  computed  from  the  pumping  power  and  the  flow  rate  using  the 
equation, 


Ap  ®  ^  (4,2,1) 

where,  Ap  is  the  pressure  drop  through  the  gill  (N/m^),  and 
P  and  Q  are  as  defined  before. 

4.2.1  Effect  of  Gill  Fiber  Outer  Diameter 


Figure  4.2. 1  Seawater  Flow  Rate  for  Various  Gill  Fiber  Outer  Diameters 
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Figures  4.2.1,  4.2.2,  and  4,2.3  depict  the  variation  of  the  seawater  flow  rate,  seawater 
pumping  power,  and  seawater  pressure  drop  respectively  as  a  function  of  the  gill  volume, 
The  graphs  show  this  variation  for  the  four  different  values  of  the  fiber  outer  diameter  listed 

in  Table  4.2.1.  The  calculations  were  done  with  the  other  parameters  (§,  1,  and  Cgi)  at  their 
respective  baseline  values  shown  in  Table  4,2,1.  The  fiber  outer  diameters  in  the  legends 
of  these  graphs  arc  in  units  of  meters, 

It  is  clear  from  Figure  4,2,1  that  as  the  volume  of  the  gill  is  increased  for  a  particular  value 
of  the  fiber  outer  diameter,  the  flow  rate  required  for  an  oxygen  flux  of  1  standard  liter  per 
minute  decreases,  This  is  because  a  larger  gill  packs  a  greater  surface  area  of  membrane 
and  therefore  requires  smaller  seawuter  velocities  for  the  same  flux,  It  is  also  apparent  that 
as  the  liber  outer  diameter  increases,  the  flow  rate  corresponding  to  a  particular  gill  volume 

increases.  ThiN  huppens  because  for  the  same  gill  volume  und  spacing  factor  $  the  larger 
Hirers  mean  a  smaller  membrane  surface  area,  and  hence  a  greater  required  seawater 
velocity.  It  is  Interesting  to  note  from  Figure  4,2,1,  however,  that  the  seawater  flow  rate 
reaches  an  asymptotic  minimum  for  gill  volumes  greater  than  approximately  0. 1  m  V  This 
is  reflective  of  the  fact  that  for  large  gills,  the  large  surface  area  and  high  residence  times 
combine  to  ensure  close  to  100%  oxygen  extraction.  The  asymptotic  minimum  therefore 
reflects  the  minimum  amount  of  seawater  Dow  rate  required  given  its  initial  oxygen 
concentration  and  the  oxygen  flux, 


Figure  4.2.2  Seawater  Pumping  Power  for  Various  Gill  Fiber  Outer  Diameters 

Figure  4.2,2  depicts  the  variation  of  the  seawater  pumping  power  as  a  function  of  the  gill 
volume  for  various  fiber  sizes,  As  expected,  the  pumping  power  decreases  as  the  volume 
increases.  At  the  lower  end  of  the  gill  volume  spectrum,  the  power  to  pump  seawater 
through  the  gill  goes  up  as  the  fiber  diameter  increases.  This  is  because  the  smaller  area 
presented  by  the  larger  fibers  dictates  that  larger  quantities  of  seawater  be  pumped  through 
the  gill.  At  higher  values  of  gill  volume,  however,  the  converse  is  true.  Larger  fibers  in 
this  case  lead  to  smaller  values  of  pumping  power.  This  is  because  of  the  fact  that  the 
seawater  flow  rate  has  reached  its  asymptotic  minimum  and  100%  extraction  is  being 
effected.  Therefore,  the  seawater  pumping  power  becomes  solely  dependent  on  the  space 
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between  the  gill  fibers  and  this  space  is  larger  for  larger  fibers  owing  to  the  fact  that  the 
value  for  $  is  constant. 

Figure  4.2.3  is  a  graph  showing  variation  of  the  pressure  drop  through  the  gill  as  a 
function  of  gill  volume  for  various  gill  fiber  sixes, 


Figure  4,2,3  Pressure  Drop  for  Various  Gill  Filler  Outer  Diameters 

4.2,2  Effect  of  Gill  Fiber  Spacing 

Figures  4.2.4,  4.2.3,  and  4,2.6  depict  the  variation  of  the  seawater  flow  rate,  seawater 
pumping  power,  and  the  presaure  drop  through  the  gill  respectively  as  a  function  of  the  gill 

volume  for  various  values  of  the  gill  fiber  spacing  factor  The  values  of  the  other 
variables,  namely  d0, 1,  and  C8i  were  held  constant  at  their  respective  baseline  values. 

Figure  4.2,4  depicts  the  variation  of  the  seawater  flow  rate  as  a  function  of  the  gill  volume 
for  four  different  values  of  the  gill  spacing  factor.  The  flow  rate  is  shown  to  increase  for 
larger  spacing  and  this  is  because  tne  larger  spacing  implies  smsller  membrane  surface 
areas  within  the  same  gill  volume,  Therefore  larger  seawater  velocities  have  to  be  attained 
in  order  to  facilitate  the  desired  oxygen  flux,  As  in  the  previous  case  (see  Figure  4.2.1), 
the  seawater  flow  rates  attain  an  asymptotic  minimum  when  the  gill  volume  exceeds 
approximately  0, 1  cubic  meters.  Once  again,  for  large  volumes,  the  flow  rate  is  so  low  that 
virtually  all  tne  oxygen  in  the  seawater  Is  extracted  and  therefore  the  flow  rate  is  limited 
primarily  by  the  initial  oxygen  concentration  and  the  desired  oxygen  flux. 
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Figure  4,2,4  Seawater  Flow  Rate  for  Various  Gill  Fiber  Spacing  Factors) 


Figure  4.2,5  Seawater  Pumping  Power  for  Various  Gill  Fiber  Spacing  Factors 


Figure  4.2,5  represents  the  pumping  power  variation  as  a  function  of  the  gill  volume  for 
the  various  gill  spacing  factor  values.  This  graph  is  similar  to  that  in  Figure  4.2.2,  in  that 
the  pumping  power  is  higher  for  larger  values  of  gill  spacing  when  the  gilT  volume  is  small. 
This  is  because  of  the  smaller  membrane  areas  implied  in  larger  gill  spacing  factors,  which 
in  tun*  necessitates  larger  seawater  flow  velocities  and  higher  now  rates.  At  higher  gill 
volumes,  the  converse  Ts  true,  i.e.,  the  pumping  power  is  smaller  for  larger  fiber  spacing. 
This  is  because  the  seawater  is  being  exhausted  of  all  oxygen  at  this  point  and  therefore,  the 
required  flow  rate  is  constant  regardless  of  the  spacing  factor.  The  larger  spaces  between 
the  gills  imply  reduced  drag  and  correspondingly  smaller  pumping  power  values. 
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Figure  4.2,6  Pressure  Drop  for  Various  0111  Fiber  Spacing  Factors 


Figure  4,2.6  depicts  the  variation  of  the  gill  pressure  drop  for  different  values  of  the  gill 
spacing  factor. 

4.2.3  Effect  of  (Sill  Fiber  lied  Thickness 

The  effect  of  the  gill  fiber  bed  thickness  1  is  examined  in  Figures  4.2.7,  4.2,8,  and  4,2.9 
which  represent  the  variation  of  the  seawater  flow  rate,  seawater  pumping  power  and  the 
pressure  drop  respectively  is  a  function  of  the  gilt  volume.  These  calculations  have  been 

done  with  the  other  variables,  namely  d0, 4.  »«d  C»|  at  their  baseline  values  as  shown  in 
Table  4.2.1. 


Figure  4.2.7  Seawater  Flow  Rate  for  Various  Gill  Thickness  Values 
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Figure  4.2.7  shows  the  variation  of  the  seawater  flow  rate  as  a  function  of  gill  volume  for 
various  values  of  the  gill  thickness  ranging  from  2,3  cm  to  10  cm.  The  gill  was  again  steed 
for  an  oxygen  flux  of  one  standard  liter  per  minute,  Figure  4,2,7  shows  that  the  seawater 
flow  rate  required  for  the  above  flux  decreases  with  increasing  gill  thickness,  This  can  be 
attributed  to  the  fact  that  by  increasingthe  gill  thickness,  a  greater  proportion  of  the  oxygen 
present  in  the  seawater  is  extracted,  Therefore  a  smaller  rate  of  seawater  flow  is  requited 
for  gills  with  thicker  fiber  beds,  As  in  the  previous  cases,  the  seawater  flow  rate  reaches  an 
asymptotic  value  beyond  a  certain  value  of  gill  volume.  This  asymptotic  value  again 
corresponds  to  the  situation  in  which  close  to  100%  of  the  oxygen  is  extracted  from  the 
seawater  as  described  in  the  preceding  sections, 

The  tiuxiel  for  the  oxygen  transport  across  the  membrane  assumes  a  logarithmic  profile  for 
the  concentration  difference  AC  along  the  thickness  of  the  fiber  bed,  This  means  that  the 
fillers  in  the  front  portion  of  the  gill  perform  most  of  the  work  as  far  as  oxygen  extraction  is 
concerned,  Because  of  this,  adding  to  the  thickness  of  a  gill  does  not  contribute  very  much 
to  the  extent  of  oxygen  extraction,  especially  at  higher  gill  volumes. 
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Figure  4.2.8  Seawater  Pumping  Power  for  Various  Gill  Thickness  Values 


Figure  4.2.8  depicts  the  variation  of  the  seawater  pumping  power  as  a  function  of  the  gill 
volume.  As  expected,  the  pumping  power  is  lowest  for  gills  with  the  smallest  fiber  bed 
thicknesses.  The  pumping  power  at  the  lower  extreme  of  the  volume  scale  Is,  however, 
independent  of  the  thickness,  This  is  because  the  positive  contribution  -  to  the  seawater 
pumping  power  -  of  increased  frontal  area  and  shorter  seawater  path  length  in  gills  with 
thinner  fiber  beds,  is  cancelled  by  the  larger  seawater  flow  rate  required.  However,  the 
larger  frontal  area  and  smaller  path  length  associated  with  thinner  gills  do  produce  power 
savings  in  gills  with  high  volumes.  This  is  because  the  virtual  depletion  of  the  seawater  in 
large  gills  dictates  that  the  flow  rate  be  independent  of  the  gill  thickness.  Hence,  it  can  be 
seen  from  Figure  4.2.8,  that  for  large  gills,  a  thinner  fiber  bed  would  offer  advantages  with 
regard  to  savings  in  seawater  pumping  power. 
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Figure  4.2.y  Pressure  Drop  for  Various  Gill  Thickness  Values 

4.2.4  Effect  of  Initial  Oxygen  Content  of  Seawater 

The  effect  of  the  initial  oxygen  concentration  in  seawater  is  examined  in  Figures  4.2.10 
through  4.2.12.  Concentrations  are  expressed  in  stundard  ml  of  02  per  ml  of  seuwuter. 

The  value.}  of  do.  and  1  were  held  constant  at  their  baseline  values  for  this  set  of 
calculations, 


Figure  4.2.10  Seawater  Flow  Rate  for  Various  Seawater  Oxygen  Concentrations 

Figure  4.2.10  shows  the  variation  of  the  seawater  flow  rate  required  for  one  standard  liter 
per  minute  oxygen  flux  as  a  function  of  the  gill  volume.  It  is  apparent  that  as  the  initial 
oxygen  concentration  decreases,  it  is  necessary  to  pump  larger  quantities  of  seawater 
through  the  gill.  The  asymptotic  minimum  for  the  seawater  flow  rate  is  different  for  each 
value  of  the  initial  seawater  oxygen  content.  This  is  to  be  expected  from  the  fact  that  the 
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asymptotic  minimum  represents  100%  extraction  of  oxygen  and  therefore  the  flow  rate 
necessary  to  produce  the  desired  flux  would  depend  only  on  the  initial  oxygen  content  of 
the  seuwater, 


Oltl  Design  Model  (1  «td  l/mln) 
Seawftter  Pumping  Power  vs,  0111  Volume 
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Figure  4,2, 1 1  Seawater  Pumping  Power  for  Various  Seawater  Oxygen  Concentrations 

The  variation  of  the  seawater  pumping  power  as  a  function  of  the  gill  volume  for  various 
values  of  Cg»  Is  investigated  in  Figure  4,2,11.  It  is  apparent  that  the  pumping  power 
required  for  conditions  of  lower  oxygen  concentration  is  higher  than  that  for  oxygen-rich 
environments. 


Figure  4,2,12  Pressure  Drop  for  Various  Seawater  Oxygen  Concentrations 

Figure  4.2.12  shows  the  pressure  drop  variation  as  a  function  of  the  gill  volume  for 
various  values  of  the  initial  oxygen  concentration  in  seawater. 
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4.3.0  Gill  Optimization 

The  DARPA  project  calls  for  the  utilization  of  a  gill  to  extract  oxygen  from  the  seawater  for 
use  in  a  stationary  power  source.  The  best  design  of  such  a  power  source  would  imply  the 
choice  of  a  gill  that  uses  the  least  amount  of  power  for  seawater  pumping.  Ideally  a  power 
system  would  not  require  any  seawater  pumping  at  ail,  but  would  rely  on  the  existing 
seawater  currents  for  oxygen  extraction.  However,  given  the  paucity  of  information 
regarding  the  magnitude  of  such  seawater  currents  and  their  transient  nature  it  is 
recommended  that  a  power  system  designed  around  a  gill  have  a  seawater  pump.  Since  it 
is  desired  that  this  gill  operate  at  the  lower  end  of  the  power  consumption  spectrum,  it  is 
likely  to  be  of  relatively  large  volume.  Figures  4.2.1,  4.2.5  and  4.2.8  investigated  the 

power  consumption  of  a  gill  for  various  values  of  do,  and  1.  It  was  found  from  the  first 
two  graphs  thut  larger  fibers  and  larger  interstitial  spaces  do  not  automatically  mean  higher 
pumping  powers.  Intuitively,  it  is  apparent  that  the  larger  gill  fiber  and  spacing  factor 
contribute  to  reduced  surface  area  available  for  mass  transport  in  a  given  volume,  and 
therefore,  higher  seawater  pumping  power,  This  higher  seawater  pumping  power  would 
be  associated  with  increased  seawater  flow  rates  and  in  turn  higher  seawater  velocities  to 
overcome  the  effect  of  the  reduced  surface  areas, 


However,  the  model  actually  predicts  lower  pumping  powers  for  larger  values  of  d0  and  ^ 
when  the  volume  of  the  gill  is  large.  This  anomaly  was  explained  as  being  due  to  the 
complete  extraction  of  oxygen  from  seawater.  The  exhaustion  of  the  seawater  oxygen 
content  is  due  to  the  comparatively  large  volume  of  the  gill  and  consequently  its  large 

surface  area.  However,  if  this  analysis  is  extended  for  ever  increasing  values  of  do  and 
it  becomes  apparent  that  there  should  be  a  limit  to  the  decrease  in  seawater  pumping  power. 

In  other  words,  the  power  would  decrease  as  the  values  of  do  and  £  arc  increased,  up  to  a 
certain  point,  beyond  which  the  surface  area  reduction  would  start  to  predominate  and  the 
conditions  of  100%  extraction  would  no  longer  exist. 

This  means  that  in  the  range  of  gill  volumes  considered  ideal  for  this  application,  it  is 
possible  to  optimize  the  values  of  do  and  %  so  as  to  obtain  a  design  requiring  minimum 
amount  of  power.  The  most  likely  process  of  design  would  commence  from  deciding  upon 
an  allocation  of  power  for  seawater  pumping.  Accordingly  the  gill  design  optimization 
process  is  illustrated  here  by  using  the  following  assumptions: 


Gill  Pumping  Power 
Oxygen  Flux  Desired 
Initial  Oxygen  Concentration 

Seawater  Viscosity 
Diffusivity  of  oxygen  in  seawater 
Density  of  seawater 


P  -  10W 

N  -  1  standard  1/min 

Csi  =  3  standard  ml/ml  of  seawater 

u  =  1.7621  E-06m2/s 

D  =  2.1  E-09  m2/s 

p  =  1025.8  kg/m3 


With  these  assumptions  in  mind  the  gill  design  was  optimized  using  Equations  4.1.21  and 
4,1.23.  The  process  commenced  by  solving  both  equations  simultaneously  using  the 
Newton-Raphson  method  to  find  values  for  V  and  Q  given  the  value  for  P. 
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4.3.1  Optimization  of  GiU  Volume  With  Respect  to  do 


Figure  4.3.1  Gill  Volume  as  a  function  of  do  for  fixed  Pumping  Power 


Figure  4.3.1  depicts  the  volume  of  a  gill  consuming  10  W  to  produce  1  1/min  of  oxygen 
given  the  conditions  listed  above.  The  volume  is  plotted  as  a  function  of  the  gill  fiber  outer 
diameter  for  various  values  of  the  fiber  bed  thickness  1  (units  in  the  legends  are  meters). 

The  gill  fiber  spacing  factor  £  was  kept  constant  at  a  value  of  2  for  this  analysis.  It  can  be 
seen  from  Figure  4.3.1  that  for  each  value  of  the  variable  1,  the  gill  volume  has  a  minimum 
value.  The  fiber  diameter  corresponding  to  the  minimum  gill  volume  changes  with  the  bed 
thickness.  This  graph  should  therefore  constitute  a  design  tool  for  the  selection  of  a  fiber 
size  given  the  desired  oxygen  flux  and  the  pumping  power  allocated  for  moving  the 
seawater. 

4.3.2  Optimization  of  Gill  Volume  With  Respect  to  ^ 


Figure  4.3.2  Gill  Volume  as  function  of  £  for  fixed  Pumping  Power 
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Figure  4.3.2  depicts  the  variation  of  the  gill  volume  V  as  a  function  of  the  gill  fiber  spacing 
factor  4  for  various  values  of  1  (units  in  the  legend  are  once  again  in  meters).  As  in  the  case 

with  do.  the  gill  volume  is  shown  to  have  a  minimum  value.  The  value  for  ^  at  which  this 
minimum  gill  volume  is  realized  also  depends  on  the  value  for  the  variable  1.  Together  with 
the  graph  shown  in  Figure  3.1,  this  graph  should  constitute  a  complete  gill  design  tool  with 

the  ability  to  specify  the  four  important  gill  design  parameters,  namely  V,  1,  do.  and  Jj, 
4.4.0  Conclusions 

The  design  of  a  gill  for  the  power  system  is  crucial  in  as  far  ns  maintaining  a  size  and  cost 
advantage  with  respect  to  other  technologies  is  concerned.  The  proper  design  of  a  gill 
specific  to  the  conditions  at  the  site  of  installation  should  help  keep  the  system  volume 
down  and  minimize  the  seawater  pumping  power  required.  These  goals  form  the  basis  of 
establishing  the  feasibility  of  the  system  as  well  as  comparing  its  performance 
characteristics  with  respect  to  other  technologies. 

In  designing  a  gill,  a  number  of  phenomena  have  to  be  considered  and  modelled.  It  is  also 
important  that  tne  model  for  the  gill  be  easily  implemented  and  be  self-sufficient  in  terms  of 
minimizing  the  number  of  assumptions.  The  model  proposed  in  this  study  does  just  that. 
The  twin  phenomena  of  drag  and  oxygen  diffusion  are  considered  and  implemented  into  a 
set  of  two  simultaneous  equations.  These  equations  form  a  well-posed  problem  that  can  be 
solved  without  the  need  to  make  extraneous  assumptions.  Moreover  the  equations  lend 
themselves  to  parametric  analysis  thus  helping  analyze  the  importance  of  the  various  gill 
sizing  parameters  and  also  assisting  in  the  Identification  of  suitable  values  for  these 
parameters, 

The  sensitivity  analysis  and  the  gill  optimization  exercise  undertaken  as  part  of  this  study 
have  yielded  general  purpose  spreadsheet  models  that  can  be  employed  as  part  of  larger 
system  sizing  models.  This  study  has  also  demonstrated  that  there  exists  a  region  defined 
by  values  of  the  gill  fiber  diameter  and  the  fiber  spacing  factor  that  allows  for  minimum 
volume  at  a  given  pumping  power  and  flux  rate.  This  region  is  specific  to  the  oxygen 
content  of  the  seawater  and  therefore  the  model  allows  flexibility  in  designing  gills 
specifically  for  the  in  situ  values  of  oxygen  concentration  at  the  site  of  application. 

Further  work  on  the  gill  design  is  intended  to  experimentally  verify  the  results  of  the  model 
developed  here.  Experimental  verification  shall  oe  performed  in  the  laboratory  with  the 
help  of  gill  cartridges  that  are  being  constructed  under  the  DARPA  project.  Eight  cartridges 
in  all  encompassing  two  variables  are  being  built.  Of  the  eight  cartridges,  four  shall  consist 
of  microporous  hollow  fiber  membranes,  while  four  shall  be  constricted  from  solid  hollow 
fiber  membranes.  The  solid  membrane  represents  the  coating  of  a  conventional 
microporous  membrane  with  a  proprietary  oxygen  permeable  polymer  to  help  prevent  pore 
collapse  and  liquid-liquid  contact  at  high  pressures.  Within  the  four  cartridges  of  each 

type,  two  values  of  the  fiber  spacing  factor  ^  shall  be  used.  These  values  will  be  1.5  and 
2.5.  These  variations  shall  help  verify  the  effect  of  the  fiber  spacing  factor  as  predicted  by 
the  model. 

In  the  experiments,  the  carrier  shall  be  pumped  through  the  insides  of  the  fibers  and  the 
seawater  through  the  fiber  bed  on  the  outside.  By  measuring  the  flow  rate  of  the  seawater, 
the  pressure  drop  in  the  seawater  flow  across  the  fiber  bed,  and  the  decrease  in  oxygen 
content  of  the  seawater,  the  model  shall  be  verified  and  calibrated. 
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5.0  DIRECT  FEED  FUEL  CELL  INTEGRATION  WITH  SEAWATER 
LOADER 

An  integration  of  the  solid  membrane  gill  module  with  the  direct  feed  fuel  cell  was  carried 
out  in  order  to  ascertain  that  the  two  subsystems  are  compatible  and  that  their  individual 
performance  is  not  compromised  by  their  interaction. 

The  intention  of  the  integration  was  to  replicate  the  conditions  used  in  the  previous  gill  and 
fuel  cell  experiments.  However  several  experimental  conditions  were  not  duplicated. 
First,  due  to  an  error  in  reading  a  rotometer  (stainless  steel  ball  instead  of  glass)  the  carrier 
flow  rate  through  the  gill  was  about  three  times  higher  than  had  previously  been  used. 
Second,  fuel  cell  experiments  had  not  previously  been  run  at  the  temperatures  used  in  the 
gill  experiments.  Third,  fuel  cell  experiments  had  not  previously  been  run  using  carrier 
solutions  without  supporting  electrolyte  as  used  in  the  gill  experiments.  These  differences 
are  reflected  in  the  results  of  the  integration  study. 

5.1  Fuel  Cell  Performance 

The  cell  voltage  at  4  mA/cm2  in  the  integration  experiment  was  0.19  volts  as  compared  to 
0.25  in  the  reference  experiment.  The  difference  can  be  accounted  for  in  part  by  the 
differences  in  temperature  and  conductivity.  The  difference  in  temperature  accounts  for 

approximately  15  mV  (AV  =  0.25  (  1  -  (277'K/  295*K)).  The  difference  in  conductivity 
accounts  for  about  13  mV  in  cell  iR.  The  difference  is  estimated  by  correcting  the 
measured  conductivities  to  5*C  and  42  mmho/cm  at  16*C  is  about  33  mmho/cm  at  5*C). 

Then  multiolying  the  measured  iR  by  the  inverse  ratio  of  the  conductivities  (AV  =  0.02  (1 
-  (33/90)'  i  he  difference  in  conductivity  will  however  result  in  a  larger  difference  in  cell 
voltage  thj.i  the  iR  difference  alone  since  the  reaction  distribution  in  the  felt  electrode  will 
be  different.  This  additional  effect  is  difficult  tw  estimate  but  could  very  well  account  for 
the  additional  32  mV  difference  in  cell  performance.  Subsequent  fuel  cell  experiments  will 
examine  the  effect  of  conductivity  on  cell  performance. 

5 . 2  Gill  Performance 

The  gill  performance  in  terms  of  maximum  oxygen  transport,  n  iliters  per  minute  per 
square  meter,  was  found  to  be  3.6  versus  2.6  in  a  previous  expr  ent  with  the  same  sea¬ 
water  flow  rate  (1.34  liters  per  minute  per  square  meter).  Howe  the  carrier  flow  rate  in 
the  previous  experiment  was  0.134  lpm/rn2  versus  0.323  Ipm/m2  n  he  integration.  The 
difference  in  oxygen  flux  can  perhaps  be  e/plained  by  the  difference  in  low  rate. 

The  maximum  oxygen  flux  was  measured  while  the  carrier  was  becoming  oxygenated  be¬ 
fore  steady  state  system  operation  was  achieved.  After  steady  state  operation  was  achieved 
the  oxygen  balance  calculated  on  the  basis  of  the  cell  current  and  the  measured  inlet  and 
outlet  partial  pressures  of  oxygen  in  seawater,  as  shown  in  Figure  5.1,  was  within 
experimental  error. 

5.3  Conclusions 

This  experiment  indicates  that  the  gill  and  fuel  cell  systems  can  be  successfully  integrated 
without  unexpected  results.  Further  development  on  both  systems  can,  therefore,  continue 
before  they  are  integrated  and  characterized  together  without  fear  that  such  development 
will  be  undone  by  the  integration. 
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6.0  ALUPOWER-ALWATT  HYDROGEN  GENERATOR 

6.1  Purpose 

This  experiment  was  to  check  the  power  output  and  hydrogen  output  of  the  Aiwatt 
Aluminum-Seawater  Battery  at  the  operating  temperature  of  zero  degrees  Celsius. 

6.2  Experimental  Setup 

The  Aiwatt  battery  was  suspended  in  thirty  gallons  of  Instant  Ocean  mixed  to  a  refractive 
index  corresponding  to  33  PPT.  The  seawater  tank  was  a  fifty  gallon  trash  can  with  a  glass 
wool  insulating  blanket  on  the  outside.  The  hydrogen  exit  of  the  battery  runs  to  a  Hoechst 
Celanese  20201040  polypropylene  membrane  cartridge  containing  three  square  meters  of 
surface  area  with  40  percent  porosity.  The  membrane  module  allows  the  hydrogen  gas  to 
flow  through  the  membrane  while  leaving  behind  the  foam  created  by  the  battery.  A  return 
form  drain  tube  is  used  so  the  foam  can  return  to  the  seawater  battery.  A  Serna  Allborg 
mass  flow  meter  with  a  0-200  ml/mi n  range  was  used  to  measure  the  flow  rate.  The  flow 
meter  was  calibrated  with  oxygen  gas  and  corrected  for  hydrogen  gas.  It  should  be  noted 
that  the  hydrogen  was  being  produced  at  the  temperature  of  the  seawater  tank  and  it  was 
being  measured  at  room  temperature. 

A  Hewlett  Packard  6825A  power  supply  was  used  as  the  power  sink  for  the  battery  output. 
Voltage  was  read  from  the  voltage  drop  across  the  power  sink  and  current  was  recorded 
from  a  shunt  resistor  (50  mV  equals  5  amps).  Data  for  the  gas  output,  cell  voltage  and  cell 
current  was  collected  using  the  Aquanautics  Ducksoup  data  collection  system. 

For  future  battery  use  it  should  be  noted  that  the  output  leads  to  the  battery  are  incorrectly 
colored  so  that  black  is  positive  and  red  is  negative. 

6.3  Procedure 

Four  sets  of  data  were  taken  during  this  experiment: 

•  The  battery  was  run  initially  in  a  short  circuit  configuration  to  evaluate  the  start-up 
profile  and  to  see  how  much  power  could  be  extracted  from  the  system.  The  time 
of  this  experiment  was  about  three  hours. 

•  The  second  set  of  experiments  was  a  current  vs.  voltage  profile  run  at  4.4  degrees 
Celsius. 

•  Next  the  cell  was  run  at  a  constant  current  based  upon  the  hydrogen  generation 
output  of  the  battery.  This  was  to  see  how  the  system  ran  at  the  desired  design 
point  based  on  hydrogen  output.  Here  the  voltage  profile  as  well  as  the  hydrogen 
output  was  measured  over  time.  This  experiment  was  run  for  19  hours.  The 
temperature  rose  6.2  degrees  Celsius  over  the  course  of  the  data  collection. 

•  The  last  set-up  change  was  a  repeat  of  the  second  experiment  where  the  voltage 
current  profile  is  examined.  For  the  second  current  voltage  profile  the  seawater 
temperature  had  risen  to  10.6  degrees  Celsius. 
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6.4  Results 

The  seawater  battery  was  first  run  in  a  short  circuit  configuration  with  no  voltage  and 
maximum  current.  This  phenomenon  is  shown  in  Figure  6.1. 

Temperature  seems  to  have  the  largest  effect  on  the  hydrogen  output  and  current  density  of 
the  cell.  Figure  6.2  shows  the  results  of  the  second  and  forth  experiments.  In  Figure  2, 
the  current-voltage  profile  is  shown  for  two  different  temperatures.  Here  the  cell  voltage 
for  zero  current  is  the  same.  However,  as  the  cell  approaches  a  short  circuit,  the  current 
output  from  the  battery  is  more  than  twenty-five  percent  higher. 

In  the  third  set  of  data,  the  voltage  is  recorded  at  a  constant  current  for  nineteen  hours. 
During  this  time  the  current  is  fixed  at  0.23  amps.  The  cell  voltage  starts  at  1.58  volts  and 
rises  to  1.66  volts.  The  corresponding  temperature  change  for  the  same  period  is  4.4 
degrees  Celsius  to  10.6  degrees  Celsius.  The  hydrogen  output  corresponds  to  the  current 
and  does  not  seem  to  be  temperature  related,  except  that  temperature  does  affect  the  current 
voltage  profile. 

An  example  of  this  is  depicted  in  Figure  6.3.  This  is  the  hydrogen  vs.  current  profile  for 
the  two  different  temperatures.  It  should  be  noted  that  the  time  period  of  data  collection  is 
different  for  the  two  experiments.  This  is  represented  in  the  data  scatter  of  the  lower 
temperature  experiment  in  Figure  6.3.  Here  the  data  for  the  first  experiment,  at  the  lower 
temperature,  was  averaged  over  a  shorter  time  period  than  the  last  set  of  data. 

6.4  Conclusions 

Once  the  cell  has  reached  its  initial  peek  operating  performance  for  a  given  temperature,  the 
temperature  of  the  seawater  seems  to  be  the  biggest  variable  in  the  battery  performance. 
The  initial  start  up  of  time  of  the  cell  seems  to  slow,  aoout  three  hours.  The  temperature 
change  of  the  seawater  during  the  start-up  period,  does  not  account  for  much  of  the  rise  in 
the  current  during  the  start-up  period. 
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7.0  CARRIER 

7.1  Summary 

6  EOCs  were  operational,  3  with  pure  O2  output  The  best  lifetime  is  6  weeks. 

7.2  Objectives 

1 .  Install  one  more  MP  cell-sized  EOC  and  complete  installation  of  "kinetics  EOC". 
Operate  at  least  3  MP  cell  systems  in  the  pure  O2  unload  mode  and  modify 
unloaders  to  optimize  for  this  output  mode. 

Status:  completed  except  for  kinetics  rig  for  which  there  has  been  no  demand. 

2 .  Repeat  pure  O2  lifetime  study  on  23SP.  Compare  with  an  air  unload  control. 

Status:  experiments  nearing  completion. 

3 .  Continue  to  screen  at  least  one  new  carrier  per  month. 

Status:  completed.  Five  new  compounds  screened  in  EOC  this  quarter. 

4 .  Purify  carriers  for  current  efficiency  comparison  purposes. 

5 .  Design  and  build  in  conjunction  with  L'Air  Liquide  long-term  EOCs. 

Status:  on  hold. 

6 .  Test  new  GDE  cell  system. 

Status:  on  hold. 

7.3  Other  Achievements 

Preparation  of  a  new  mediator  has  been  accomplished  and  the  compound  functions  well.  It 
has  shown  good  capability  for  catalysis  with  23SP  and  has  enabled  the  power  requirement 
to  be  reduced  to  less  than  150  WminL'1. 

1 .  We  now  have  six  operating  EOCs,  two  of  which  have  25  cm^  cells  and  four  have 
100  cm2  MP  cells.  Of  the  four  larger  size  systems,  three  are  capable  of  continuous 
operation  in  the  pure  O2  unload  mode  with  reasonable  accuracy.  We  have  two  10 
ml/min  and  one  20  ml/min  flowmeters  and  two  further  50  ml/mi n  flowmeters. 
Since  the  accuracy  is  not  good  below  10%  full  scale,  our  lower  limit  for  output  is 
thus  1  ml/min.  The  present  unloaders  appear  to  operate  satisfactorily  at  the  current 
oxygen  outputs.  Further  design  changes  will  be  implanted. 

2 .  One  experiment  has  been  completed  using  23SP  plus  mediator  in  an  EOC  in  the  air 
unload  mode.  The  carrier  lasted  for  six  weeks  before  failing.  Another  experiment 
has  been  underway  for  four  weeks  using  similar  carrier  conditions  but  with  a  pure 
O2  unloading  mode.  The  experiment  is  not  yet  complete. 

3 .  During  the  period,  3  new  compounds  have  been  screened,  2SP,  3SP  and  232NP. 

4 .  Experiments  with  23SP  have  been  performed  where  the  HQ  salt  was  recrystallized 
four  times  and  carefully  dried.  This  material  gave  continuous  performance  with  an 
electron  count  of  less  than  3.  Another  experiment  using  free  base  23SP  gave  an 
electron  count  of  closer  to  4.  This  results  in  poorer  power  requirements  and  poorer 
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purity  when  the  EOC  is  run  in  the  pure  O2  unload  mode.  The  steady-state  electron 
count  lifetime  does  not  appear  to  be  affected,  however. 

A  new  mediator  has  been  screened.  It  performs  as  expected  with  33SP  giving  the 
same  performance  as  earlier  mediators.  When  used  with  23SP  the  power 
requirements  dropped  to  less  than  150  watt.min/L  for  almost  four  weeks  and  the 
useful  lifetime  of  23SP  was  extended  to  6  weeks  from  4  weeks. 

During  the  course  of  support  work  for  fuel  cell  experiments,  some  experiments 
with  anions  other  than  chloride  were  tried.  The  use  of  fluoride  seemed  to  be 
acceptable,  although  the  conductivity  was  lower.  The  use  of  perchlorate  was 
disastrous,  the  solution  rapidly  polymerizing  to  give  insoluble  material.  Both 
carbonate  and  acetate  were  tried.  The  latter  functioned  in  the  fuel  cell  but  in  the 
EOC  a  voltage  of  1.2  volts  was  required  to  obtain  oxygen  from  23SP  as  opposed  to 
0.8  -  0.9  V  with  chloride.  This  indicates  binding  by  acetate  to  the  cobalt  (III)  form. 
Such  observations  are  in  line  with  our  previous  observations  with  ferrocene 
carboxylic  acids  and  phosphate. 

7.4  Objectives  and  Workplan  for  Next  Quarter 

1 .  Design,  build  and  operate  a  small  volume  EOC  requiring  less  than  50  ml.  of 
solution.  Design  and  build  new  loaders  and  unloaders  to  reduce  present  volumes 
by  50%. 

2 .  Design,  build  and  operate  a  battery  of  fuel  cell  test  apparati  to  evaluate  carriers  for 
fuel  cell  operation. 

3 .  Screen  at  least  one  new  compound  per  month  in  the  EOC  and  one  per  month  in  fuel 
cells.  Operate  at  least  three  in  the  pure  O2  unload  mode  and  achieve  the  longevity 
goal  of  6  weeks  operation  at  >90%  oxygen  purity  at  less  than  200  watts  per  minute 
per  liter. 

4.  Evaluate  different  electrode  treatment  procedures  with  respect  to  electrochemical 
kinetics.  Prepare  electrodes  modified  with  catalysts  and  test.  Initiate  kinetics 
program  with  C.  Amatore,  obtain  and  operate  basic  simulation  software. 

5.  Continue  evaluation  of  carriers  by  cyclic  voltammetry  for  stability  and  fuel  cell 
function. 

6.  Revamp  reporting  procedures  for  EOC  operations  and  carrier  evaluations.  Ensure 
that  all  EOC  databases  are  up  to  date  and  complete. 
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